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Nomenclature 


2D: 2 dimensional 

4D: 4 dimensional 

ADS-B: Automatic Dependence Surveillance Broadcast 

BOD: Bottom-Of-Descent 

CAS: Calibrated Airspeed 

DTG: Distance-To-Go 

MSL: Mean Sea Level 

RF: Radius-to-Fix 

STAR: Standard Terminal Arrivals 

TAS: True Airspeed 

TCP: Trajectory Change Point 

TOD: Top-Of-Descent 

TTG: Time-To-Go 

VTCP: Vertical Trajectory Change Point 

Subscripts 

Subscripts associated with waypoints and TCPs, e.g., TCP 2 , denote the location of the waypoint or TCP in 
the TCP list. Larger numbers denote locations closer to the end of the list, with the end of the list being 
the runway threshold. Subscripts in variables indicate that the variable is associated with the TCP with 
that subscript, e.g., Altitude 2 is the altitude value associated with TCP 2 . 
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Units and Dimensions 

Unless specifically defined otherwise, units (dimensions) are as follows: 

time: seconds 

position: degrees, + north and + east 

altitude: feet, above MSL 

distance: nautical miles 

speed: knots 

track: degrees, true, beginning at north, positive clockwise 
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Abstract 


This document describes an algorithm for the generation of a four 
dimensional trajectory. Input data for this algorithm are similar to an 
augmented Standard Terminal Arrival (STAR) with the augmentation in 
the form of altitude or speed crossing restrictions at waypoints on the 
route. This version of the algorithm accommodates constant radius turns 
and cruise altitude waypoints with calibrated airspeed, versus Mach, 
constraints. The algorithm calculates the altitude, speed, along path 
distance, and along path time for each waypoint. Wind data at each of 
these waypoints are also used for the calculation of ground speed and 
turn radius. 

Introduction 

Concepts for self-spacing of aircraft operating into airport terminal areas have been under development 
since the 1970's (refs. 1-20). Interest in these concepts has recently been renewed due to a combination of 
emerging, enabling technology (Automatic Dependent Surveillance Broadcast data link, ADS-B) and the 
continued growth in air traffic with the ever increasing demand on airport (and runway) throughput. 
Terminal area self-spacing has the potential to provide an increase in runway capacity through an increase 
in the accuracy of runway threshold crossing times, which can lead to a decrease of the variability of the 
runway threshold crossing times. Current concepts use a trajectory based technique that allows for the 
extension of self-spacing capabilities beyond the terminal area to a point prior to the top of the en route 
descent. 

The overall NASA Langley concept for a trajectory-based solution for en route and terminal area self- 
spacing is fairly simple and was originally documented in reference 21. By assuming a 4D trajectory for 
an aircraft and knowing that aircraft’s position, it is possible to determine where that aircraft is on its 
trajectory. Knowing the position on the trajectory, the aircraft’s estimated time-to-go (TTG) to a point, in 
this case the runway threshold, is known. To apply this to a self-spacing concept, a TTG is calculated for 
a leading aircraft and for the ownship. Note that the trajectories do not need to be the same. The nominal 
spacing time and spacing error can then be computed as: 

nominal spacing time = planned spacing time interval + traffic TTG. 

spacing error = ownship TTG - nominal spacing time. 

The foundation of this spacing concept is the ability to generate a 4D trajectory. The algorithm 
presented in this paper uses as input a simple, augmented 2D path definition (i.e., a traditional STAR, 
with relevant speed and altitude crossing constraints) along with a forecast wind speed profile for each 
waypoint. The algorithm then computes a full 4D trajectory defined by a series of trajectory change points 
(TCPs). The input speed (Mach or CAS) or altitude crossing constraint includes the deceleration rate or 
vertical angle value required to meet the constraint. The TCPs are computed such that speed values, Mach 
or CAS, and altitudes change linearly between them. TCPs also define the beginning and ending segments 
of turns, with the midpoint defined as a fly-by waypoint. The algorithm also uses the waypoint forecast 
wind speed profile in a linear interpolation to calculate the wind speed at the altitude the computed 
trajectory crosses the waypoint. Wind speed values are then used to calculate the ground speeds along the 
path. 

The major complexity in computing a 4D trajectory involves the interrelationship of ground speed with 
the path distance around turns. In a turn, the length of the estimated ground path and the associated turn 
radius will interact with the waypoint winds and with any change in the specified speed during the turn, 
i.e., a speed crossing-restriction at the waypoint. Either of these conditions will cause a change in the 
estimated turn radius. The change in the turn radius will affect the length of the ground path which can 

1 



then interact with the distance to the deceleration point, which thereby affects the turn radius calculation. 
To accommodate these interactions, the algorithm uses a multi-pass technique in generating the 4D path, 
with the ground path estimation from the previous calculation used as the starting condition for the 
current calculation. 

Algorithm Overview 

The basic functions for this trajectory algorithm are shown in figure 1. Figure 1 also contains logic and 
some simple calculations that are not included in the body of this document. Also note that waypoints are 
considered to be TCPs but not all TCPs are waypoints. 

For the 2D input, the first and last waypoints must be fully constrained, i.e., have both a speed and 
altitude constraint defined. With the exception of the first waypoint, which is the waypoint farthest from 
the runway threshold, constraints must also include a variable that defines the means for meeting that 
constraint. For altitude constraints, this is the inertial descent angle; for speed constraints, it is the air 
mass CAS deceleration rate. A separate, single Mach-to-CAS transition speed (CAS) value may also be 
input for profiles that involve a constant Mach / CAS descent segment. Additionally, an altitude / CAS 
restriction (e.g., in the U.S., the 10,000 ft / 250 kt restriction) may also be entered. 

The algorithm computes the altitude and speed for each waypoint. It also calculates every point along 
the path where an altitude or speed transition occurs. These points are considered vertical TCPs (VTCPs). 
TCPs also define the beginning and ending segments of turns, with the midpoint defined as a fly-by 
waypoint. Turn data are generated by dividing the turn into two parts (from the beginning of the turn to 
the midpoint and from the midpoint to the end of the turn) to provide better ground speed (and resulting 
turn radius) data relative to a single segment estimation. A fixed, average bank angle value is used in the 
turn radius calculation. The algorithm also uses the forecast wind speed profile for a waypoint in a linear 
interpolation to calculate the wind speed at the altitude the computed trajectory crosses the waypoint (if 
the crossing altitude is not at a forecast altitude). For non-waypoint TCPs, the generator uses the forecast 
wind speed profile from the two waypoints on either side of the TCP in a double linear interpolation 
based on altitude and distance (to each waypoint). Of significant importance for the use of the data 
generated by this algorithm is that altitude and speeds (Mach or CAS) change linearly between the TCPs, 
thus allowing later calculations of DTG or TTG for any point on the path to be easily performed. 
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Trajectory calculation: 

2D input data, crossing data, 
and wind forecast data / 


Preprocess RF legs : For each waypoint identified as an RF turn initiation point, generate 
one or two pseudo fly-by waypoints for the center of the turn. 


Save the crossing angles, descent speeds, and the initial Mach: Save the altitude 


crossing angles for each waypoint, the Mach for the first waypoint, the descent Mach, 
and descent CAS values. These values may be overwritten and will need to be reset to 
their original values for each iterative loop. 


Generate the initial tracks and distances : Using great-circle calculations, determine the 
distances and ground tracks between waypoints. Calculate the DTG for each waypoint. 


Initialize the waypoint turn data : Waypoints that have more than a 3 degree change in 
ground track from the previous waypoint are considered tum-waypoints. Mark each as 
a tum-waypoint and insert a turn-entry and turn-exit TCP on each side of this waypoint. 



Reinstate the descent speeds : Restore the descent Mach and descent CAS values to their 
original values. 


Compute the TCP altitudes : Beginning at the runway (the last waypoint) work 
backwards and compute the altitude at each prior TCP. If an altitude is computed to be 
reached prior to the previous TCP, insert a new altitude TCP. 


Copy crossing angles : Beginning at the runway, for TCPs that do not have crossing 
angles, copy the downstream angle into this TCP. 


Meet cruise CAS waypoint restriction: If required, change the descent Mach if there is a 
cruise CAS restricted waypoint and the computed speed is above the required crossing 
speed. 


Add descent Mach waypoint: If the descent is to be performed initially in a Mach mode 
and the descent Mach is different than the cruise Mach, change the descent waypoint 
Mach as necessary and add any required, additional TCPs. 


Compute the Mach-to-CAS TCP: If required, compute the Mach-to-CAS transition 
altitude. Compute the DTG to this altitude and insert the Mach-to-CAS TCP. 


continued 
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Figure 1 (continued). Basic functions. 
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Figure 1 (continued). Basic functions. 
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Algorithm Input Data 

The algorithm takes as input a list of waypoints, their trajectory-specific data, and associated wind 
profile data. The list order must begin with the first waypoint on the trajectory and end with the runway 
threshold waypoint. The trajectory-specific data includes: the waypoint's name and latitude / longitude 
data, e.g., Latitude 2 and Longitude 2 , an altitude crossing restriction, if one exists, and its associated 
crossing angle, e.g., Crossing Altitude 2 and Crossing Angler, and a speed crossing restriction (Mach or 
CAS), if one exists, and its associated CAS rate, e.g., Crossing CAS 2 and Crossing Rate 2 . A value of 0 as 
an input for an altitude or speed crossing constraint denotes that there is no constraint at this point. A 
Crossing Mach may not occur after any non-zero Crossing CAS input. The units for Crossing Rate are 
knots per second. 

In this algorithm, a radius-to-fix (RF) segment is indicated by the addition of a center-of-tum position, 
e.g., Center of Turn Latitude 2 and Center of Turn Longitude 2 , for the input waypoint at the initiation of 
the turn. Additional requirements for the RF segment are provided in a subsequent section. 

To accommodate a descent from the cruise altitude, a Mach value, Mach Descent Mach, may be 
specified that is different from the cruise Mach value. A CAS value may also be specified for the Mach- 
to-CAS transition speed, Mach Transition CAS, during the descent. Additionally, a CAS speed limit at a 
defined altitude may also be included. In the U.S., this would typically be set to 250 kt at 10,000 ft. 

For the wind forecast, a minimum of two altitude reports (altitude, wind speed, and wind direction) 
should be provided at each waypoint. The altitudes should span the estimated altitude crossing at the 
associated waypoint. The algorithm assumes that the input data are valid. 

Internal Algorithm Variables 

The significant variables computed by this algorithm are: 

Altitude the computed altitude at the TCP 

CAS the computed CAS at the TCP 

DTG the computed, cumulative distance from the runway 

Ground Speed the computed ground speed at the TCP 

Ground Track the computed ground track at the TCP 

Mach the computed Mach at the TCP 

TTG the computed, cumulative time from the runway 

Additionally, the algorithm denotes TCPs in accordance with how they are generated. TCPs are identified 
as: 

• Input, from the input waypoint data; 

• An internally generated, radius-to fix (RF) center of turn waypoint; 

• Turn-entry, identifying a TCP that marks the start of a turn; 
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Turn-exit, identifying a TCP that marks the end of a turn; and 




• Vertical TCPs (VTCPs), denoting a change in the altitude or speed profile. 

TCPs may also be marked with a vertical identifier denoting one of the following: 

• Altitude, denoting a change in the descent angle; 

• Speed, denoting a change in the CAS or Mach; 

• Top of descent point, TOD; 

• Altitude CAS restriction, denoting a speed change due to a speed restriction at a specific altitude, e.g., 
250 kt at 10,000'; and 

• Mach-to-CAS, denoting the Mach-to-CAS transition point. 

TCPs are also denoted relative to the associated primary speed value, i.e., the crossing speed is Mach or 
CAS derived. 

There are also several input variables that may become overwritten within the algorithm that are 
required to be restored for subsequent calculation cycles within the algorithm. These variables include the 
following: 

• Saved Mach Descent Mach, which is the saved input value of Mach Descent Mach. 

• Saved Mach Transition CAS, which is the saved input value of Mach Transition CAS. 

• Saved Mach at First Waypoint, which is the saved input Mach value for the first waypoint, i.e., 
Crossing Mach First wavpoim, assuming that one exists. 
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Description of Major Functions 

The functions shown in figure 1 are described in detail in this section. The functions are presented in 
the order as shown in figure 1. Secondary functions are described in a subsequent section. In these 
descriptions, the waypoints, which are from the input data and are fixed geographic points, are considered 
to be TCPs but not all TCPs are waypoints. Nesting levels in the pseudo-code description are denoted by 
the level of indentation of the document formatting. Additionally, long sections of logic may end with end 
of statements to enhance the legibility of the text. 

Preprocess RF Legs 

A radius-to-fix (RF) turn segment is a constant radius turn between two waypoints, with lines tangent 
to the arc around a center of turn point (fig. 2). This function determines if a valid RF turn exists, and if 
so, calculates a pseudo-waypoint relative to the center-of-tum point and inserts it into the waypoint list. 
The calculated pseudo-waypoint then allows the remainder of the turn calculations performed by this 
algorithm to be processed as a standard turn. This function is performed in the following manner: 


TCP 



TCP l (RF turn initiation waypoint) 
Pseudo -waypoint 


TCP t + j (turn exit) 


Center of turn point 


TCP, 


i+2 


Figure 2. Example of an RF turn. 


error = false 

Big Turn Error = false 

A set of RF turn waypoints is identified by the inclusion of a non-zero value for the latitude and 
longitude for the center of turn point in the data for the RF turn initiation waypoint. Because three 
waypoints are needed in an RF turn calculation, two each for the determination of the inbound and 
outbound track angles, testing is only performed to the number of the last waypoint - 2. 

for (i = index number of the first waypoint + 1; i < index number of the last waypoint - 2; i = i + 1) 

Determine if this is an RF turn waypoint via the inclusion of the turn center's latitude and 
longitude data. 

if ((Center Of Turn Latitude, f 0) and (Center Of Turn Longitude , f 0 )) 

Determine the turn direction. 

ai = arctangent2(sine( Longitude, - Longitude,^) * cosine(Latitudej), cosinef Latitude,^) * 
sine(Latitudei) - sine(Latitude h t ) * cosine(Latitudei) * cosine(Longitudej - Longitude 



a 3 = arctangent2( sine(Longitude i+ i - Longitude) * cosine(Latitude i+ i), cos in e( Latitude,) * 
sine(Latitude i+ i) - sine(Latitudei) * cosine (Latitude i+ i) * cosine(Longitude i+ i - 
Longitude,)) 

deltax = DeltaAngleta /, a 3 ) 

where the secondary function DeltaAngle is described in a subsequent section. 

If deltax is positive, this is a right-hand turn. 
if (deltax > 0) TurnSign = 1 
else TurnSign = -1 

Calculate the instantaneous angle at the ending waypoint. 

a 2 = arctangent2( sine(Longitude i+ i - Center Of Turn Longitude ) * cosine(Latitude i+1 ), 
cosine(Center Of Turn Latitude ) * sine(Latitude i+1 ) - sine(Center Of Turn Latitude) * 
cosine(Latitude i+ i) * cosinet Longitude, - 1 - Center Of Turn Longitude)) + 
TurnSign * 90 

Adjust a 2 such that 0 > a 2 > 360 
deltaa = DeltaAngleta h a 2 ) 

Correct the deltaa value if it is in the wrong direction. 
if ((TurnSign > 0) and (deltaa < 0)) 
deltaa = deltaa + 360 
else if ((TurnSign < 0) and (deltaa > 0)) 
deltaa = deltaa - 360 

If the turn is greater than 170°, break it into two parts so that the standard turn calculations 
can be performed. 

if (\deltaa\ > 170) BigTurn = true 

If the turn is less than 3° or more than 260°, it is in error. 

if ((\deltaa\ < 3) or (\deltaa\ > 260)) error = true 

Perform a center-of-tum test. 

if (error = false) 

The radius for point 1 must equal the radius for point 2. 

/•/ = arccosine( sine(Center Of Turn Latitude) * sine (Latitude) + cosine(Center 
Of Latitude) * cosine(Latitude) * cosine(Center Of Turn Longitude , - 
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Longitude ,) ) 

r 2 = arccosine( sine(Center Of Turn Latitude ) * sine (Latitude^ i) + 
cosine(Center Of Turn Latitude j * cosine(Latitude i+ i) * 
cosine(Center Of Turn Longitude , - Longitude i+ i) ) 

The radii are considered not equal if the difference is greater than 200 ft. The overall RF 
leg is considered in error if the turn radius is greater than 10 nmi. 

if ((] ri - r 2 1 > (200 / 6076)) or (rj > 10)) error = True 

if (error = false) 

If the turn is greater than 170°, generate two waypoints, otherwise, just generate one 
waypoint. 

if (BigTurn) n = 2 

else n = 1 

a = TurnSign * 90 

for (k = 1; k<n; k = k + 1) 

Calculate the pseudo-RF waypoint. 

The following is the angle from the turn center toward the pseudo waypoint. 
a 3 = aj - a 

Adjust a 2 such that 0>a3> 360 
if (BigTurn) 

if (k = 1) a Ib = a 3 + 0.25 * deltaa 
else a ib = as + 0.75 * deltaa 

else 

Just one new waypoint, split the turn in half. 
a ib = 03 + 0.5 * deltaa 
Adjust a ib such that 0>ai b > 360 

If (k = 1) 

RadialRadiallntercept (Latitude j, Longitude h a/, 

Center Of Turn Latitude „ Center Of Turn Longitude „ ai h , 

Latitude,-), Longitude,.)), 
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noting that Latitude ^ and Longitude,-/ arc returned values. 


else 


RadialRadialIntercept(Latitude i+ i, Longitude i+I , a 2 + 180, 

Center Of Turn Latitude Center Of Turn Longitude a lb , 

Latitude r f, Longitudef , 

The new waypoint is inserted at location i+1 in the waypoint list. This inserted 
waypoint will appear as an input waypoint to the remainder of the algorithm. The 
waypoint is inserted between waypoint, and waypoint, +/ from the original list. The 
function InsertWaypoint should be appropriate for the actual data structure 
implementation of this function. 

InsertWaypoint(i + 1) 

Note that Wpt i+ i is the newly created waypoint. 

Mark Wpt i+ i as though it was an input waypoint and give it a unique name. 

Also marking this waypoint as a special, RF turn center waypoint. This special 
marking is used in subsequent sections to denoted that the center-of-tum point has 
already been calculated. 

Wpti+i = rf-turn-center 

Latitude i+ i = Latitude,/ 

Longitude i+I = Longitude ,.f 

Copy the wind data from Wpt i} the RF initiation waypoint, to Wpt i+ i, the pseudo- 
waypoint. 

Save the center of turn data. The Turn Data values are associated with each waypoint 
or TCP record and contain, if appropriate, data relating to turn conditions for that 
TCP. 

Turn Data Center Latitude i+ i = Center Of Turn Latitude , 

Turn Data Center Longitude i+ i = Center Of Turn Longitude , 

Increment i because a waypoint was added and the new waypoint at i + / should not 
be processed again. 

i = i + 1 

end of for (k = 1; k<n; k = k + 1) 
end of if (error = false) 

end of if ((Center Of Turn Latitudei f 0) and (Center Of Turn Longitude , f 0 )) 
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end of for (i = index number of the first waypoint + 1; ...) 

Generate Initial Tracks and Distances 

This is an initialization function that initializes the Mach Segment flag, denoting that the speed in this 
segment is based on Mach, and calculates the point-to-point distances and ground tracks between input 
waypoints. Great circle equations are used for these calculations, noting that the various dimensional 
conversions, e.g., degrees to radians, are not shown in the following text. 

Generate the initial distances, the center-to-center distances, and ground tracks between input 
waypoints 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 

Start with setting the Mach segments flags to false. 

Mach Segment , = false 

Compute the waypoint-center to waypoint-center distances. 

if (i = index number of the first waypoint) Center to Center Distance, = 0 

else 


Center to Center Distance , = 

arccosine( sine(Latitudej.i) * sine(Latitudej) + cosine(Latitudei-i) * cosine(Latitude ) * 
cosine(Longitudei-i - Longitude ,) ) 

Ground Tracks) = 

arctangent2( sine(Longitudej - Longitude,- f) * cosine(Latitudei), cosine(Latitudei-i) * 
sine(Latitudej) - sine(Latitudei-i) * cosine(Latitudei) * cosine(Longitudej - 
Longitude i.])) 

end of for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 
Now set the runway's ground track. 

Ground Tr acki as t waypoint Gr ound Tr acki as t waypoint - i 
The cumulative distance, DTG, is computed as follows: 

DTG, ast waypoint 0 

for (i = index number of the last waypoint; i > index number of the first waypoint; i = i - 1) 

DTGpi = DTGi + Center to Center Distance, 
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Initialize Waypoint Turn Data 

The Initialize Waypoint Turn Data function is used to determine if a turn exists at a waypoint and if so, 
inserts turn-entry and turn-exit TCPs. Waypoints that have more than a 3 degree change in ground track 
between the previous waypoint and the next waypoint are considered tum-waypoints. The function is 
performed in the following manner: 

i = index number of the first waypoint + 1 

Last Track = Ground Trackf irs , way pomt 

Note that the first and last waypoints cannot be turns. 

while (i < index number of the last waypoint) 

Track Angle After = Ground Track , 
a = DeltaAngle(Last Track, Track Angle After) 

Check for a turn that is greater than 170 degrees. 
if(\a\ > 170) 

Set an error and ignore the turn. 

Mark this as an error condition, 
a = 0 

If the turn is more than 3 -degrees, compute the turn data. 

if (l«l > V 

half turn = a / 2 

Track Angle Center = Last Track + half turn 

This is the center of the turn, e.g., the original input waypoint. 

Ground Track , = Track Angle Center 
Turn Data Trackf = Last Track 
Turn Data Track2j = Track Angle After 

If this is not an RF turn, then the turn radius needs to be calculated. 
if (Wpti f rf-turn-center) Turn Data Turn Radiusi = 0 
Turn Data Path Distance , = 0 
Insert a new TCP at the end of the turn. 
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The new TCP is inserted at location i+1 in the TCP list. The TCP is inserted between TCP, 
and TCP,-/ from the original list. The function Insert Way point should be appropriate for the 
actual data structure implementation of this function. 

InsertWaypoint(i + 1) 

Note that TCP i+1 is the new TCP. 

TCP i+i = turn-exit 

DTG i+I = DTG i 

Ground Track i+I = Track Angle After 
The start of the turn TCP is as follows, 

InsertWaypoint( i) 

TCP i = turn-entry 

Note that the original TCP is now at index i+1. 

DTG, = DTGi+i 
Ground Track , = Last Track 
Last Track = Track Angle After 
i = i + 2 

end of if (\a\ > 3) 
else Last Track = Ground Track i 
i = i + 1 

end of while (i < index number of the last waypoint) 

Effectively, this function: 

- Marks each tum-waypoint and sets its ground track angle to the computed angle at the midpoint of 
the turn. 

- Inserts a co-distance turn-entry TCP before this tum-waypoint with the ground track angle for this 
turn-entry TCP set to the value of the inbound ground track angle. 

- Inserts a co-distance turn-exit TCP after this tum-waypoint with the ground track angle for this 
turn-exit TCP set to the value of the outbound ground track angle. 

An example illustrating the inserted turn-start and tum-end TCPs is shown in figure 3. 
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Turn waypoint, Ground Trackj = 105° 

v- Turn-exit, Ground Track: + j = 120 ° 


Ground Track j_? = 90° 

Turn-entry, Ground Track j_ ] = 90° 
DTGj_j = DTG j+ j = DTGj 


Ground Trackj + 2 = 120 ° 
Figure 3. Initialized turn waypoint. 


Compute TCP Altitudes 

Beginning with the last waypoint, the Compute TCP Altitudes function computes the altitudes at each 
previous TCP and inserts any additional altitude TCPs that may be required to denote a change in the 
altitude profile. The function uses the current altitude constraint ( TCP , in fig. 4), searches backward for 
the previous constraint (TCP/.s in fig. 4), and then computes the distance required to meet this previous 
constraint. The altitudes for all of the TCPs within this distance are computed and added to the data for 
the TCPs. If the along-path distance to meet the previous constraint is not at a TCP, a new altitude VTCP 
is inserted at this distance. An example of this is shown in figure 5. In addition, if the Crossing Angle for 
a waypoint is set to -99, this denotes that the algorithm is to internally compute the Crossing Angle 
between this and the next higher, altitude constrained waypoint, noting that this option should only be 
used in situations where the relevant waypoint pairs are known to procedurally have a fixed angle 
between them. This function is performed in the following steps: 

Crossing Altitude^ = 11000ft 

Crossing Altitude l = 5300 ft 
Crossing Angle t =2.3° 


tcp i-3 tcp i-2 tcp i- 1 TCPl 


Figure 4. Input altitude crossing constraints. 


Crossing Altitude i _3 = 11000ft 


Altitude i _2 — 9560 ft 


- - ^Addude^i 7940 ft Crossing Altitude l = 55 

Crossing Angle t = 2. 

"" 


= 5300 ft 
= 2.3° 


TCP u3 New TCP 


TCP 


i-2 


TCP , 


i-1 


TCP, 


Figure 5. Computed altitude profile with TCP added. 
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Set the current constraint index number, cc, equal to the index number of the last waypoint, 

cc = index number of the last waypoint 

Set the altitude of this waypoint to its crossing altitude, 

Altitude cc = Crossing Altitude cc 

Set a flag denoting that the TOD point has not been identified 
Have TOD = false 

While (cc > index number of the first waypoint) 

If this is the TOD, mark this point. 

if Have TOD is false and Altitude cc is equal or greater than Altitude i 
Have TOD = true 
mark this as the TOD point. 

Determine if the previous constraint cannot be met. 

If (Altitude C c > Crossing Altitude C c) 

The constraint has not been made. 

If this is the last pass through the algorithm, mark this as an error condition. 

Altitude C c = Crossing Altitude cc 

Find the prior waypoint index number pc that has an altitude constraint, e.g., a crossing altitude 
(Crossing Altitude pc T 0). This may not always be the previous (i.e., cc - 1 ) waypoint. 

Initial condition is the previous TCP. 

pc = cc - 1 

while ( (pc > index number of the first waypoint) and ( (TCP pc f input waypoint) or 
(Crossing Altitude pc = 0)) ) pc = pc - 1 

Save the previous crossing altitude, 

Prior Altitude = Crossing Altitude pc 

Save the current crossing altitude (Test Altitude) at TCP,, and the descent angle (Test Angle) 
noting that the first and last waypoints always have altitude constraints and except for the first 
waypoint, all constrained altitude points must have descent angles. 

Test Altitude = Crossing Altitude cc 
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Test Angle = Crossing Angle cc 

If the Test Angle value, i.e., AUTO DESCENT ANGLE, denotes that this is angle is to be 
computed internally as a linear descent between the two altitude constrained waypoints then the 
following calculations are performed: 

if (Test Angle = A UTO DESCENT ANGLE) 

dx = DTGp C - DTG CC 

dy = Prior Altitude - Test Altitude 

Test Angle = arctangent2 (dy, 6076 * dx) 

Crossing Angle cc = Test Angle 

Test for an extreme angle, e.g., 7.5°. 

if (Test Angle > maximum allowable descent angle) mark this as an error condition. 

Compute all of the TCP altitudes between the current TCP and the previous crossing waypoint. 
k = cc 
while k> pc 

If the previous altitude has already been reached, set the remaining TCP altitudes to the 
previous altitude. 

if (Prior Altitude < Test Altitude) 

for (k = k - 1; k> pc; k = k - 1) Altitude k = Test Altitude 

Set the altitude at the last test point. 

Altitude pc = Test Altitude 

else 

Compute the distance from TCP/, to the Prior Altitude using the altitude difference 
between the Test Altitude and the Prior Altitude with the Test Angle. If there is no point at 
this distance, add a TCP at that distance. 

Compute the distance dx to make the altitude. 

dx = (Prior Altitude - Test Altitude) / (6076 * tangent(Test Angle)) 

Compute the altitude z at the previous TCP. 

z = ((DTGk-i - DTGi ) * 6076) * tangent(Test Angle) + Test Altitude 
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If there is a TCP prior to this distance or if z is very close to the Prior Altitude , compute 
and insert its altitude. 

if ( ( DTG k . i < (DTG k + dx)) or ( \z - Prior Altitude\ < some small value) ) 
if (\z - Prior Altitude\ < some small value) Altitude k -i = Prior Altitude 
else Altitude k-i = z 

Check to see if the constraint has been reached with a 100 ft tolerance; if not, set an 
error condition. 

if ( (k-1) =pc) 

if ( \Altitude pc - Crossing Altitude pc \ > 100ft ) mark this as an error condition 
Always set the crossing exactly to the crossing value. 

Altitude pc = Crossing Altitude pc 
Update the Test Altitude. 

Test Altitude = Altitude k _i 

Decrement the counter to set it to the prior TCP. 

k = k- 1 

end of if ( (DTG k _i < (DTG k + dx) ) or ( \z - Prior Altitude\ < some small value) ) 
else 

The altitude constraint is reached prior to the TCP, a new VTCP will need to be 
inserted at that point. The distance to the new TCP is, 

d = DTG k + dx 

Compute the ground track at distance d along the trajectory and save it as Saved 
Ground Track. 

Saved Ground Track = GetTrajGndTrk(d) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted between 
TCP/;./ and TCP/ from the original list. The function InsertWaypoint should be 
appropriate for the actual data structure implementation of this function. 

InsertWaypoint(k) 

Update the data for the new VTCP which is now TCP k . 
if (VSegType k = no type) VSegType k = ALTITUDE 
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DTG k — d 


Altitudek = Prior Altitude 

Add the ground track data which must be computed if the new VTCP occurs within a 
turn. The functions WptlnTurn and ComputeGndTrk are described in subsequent 
sections. 

if (WptlnTurn(k)) Ground Tracks = ComputeGndTrk(k, d) 
else Ground Track k = Saved Ground Track 

Compute and add the wind data at distance d along the path to the data of TCP k . 
GenerateWptWindProfile(d, TCP A 
Test Altitude = Prior Altitude 

Since TCP k , has now been added prior to pc, the current constraint counter cc needs 
to be incremented by 1 to maintain its correct position in the list. 

cc = cc + 1 

The function loops back to while k > pc. 

Now go to the next altitude change segment on the profile. 
cc = k 

The function loops back to while cc > index number of the first waypoint. 

Copy Crossing Angles 

The Copy Crossing Angles is a simple function that starts with the next to last TCP and copies the 
subsequent crossing angle if the current TCP does not have a crossing angle. E.g., 

for (i = index number of the last waypoint - 1; i > index number of the first waypoint; i = i - 1) 

if (Crossing Angle, = 0) Crossing Angle t = Crossing Angle i+ i 

Meet Cruise CAS Waypoint Restriction 

The Meet Cruise CAS Waypoint Restriction function changes, if required, the descent Mach if there is a 
high altitude, CAS restricted waypoint and the computed speed is above the required crossing speed for 
that CAS waypoint. 

The calling function provides as input and retains the subsequent outputs for the following variables: 
Todld, TodMach, TodMachRate, and MachCasAtTod. The variable Todld is the name of the top-of- 
descent waypoint (TOD) and is initialized as an empty string by the calling program. This Meet Cruise 
CAS Waypoint Restriction function may modify the Mach and speed change rate that occurs at the TOD, 
TodMach and TodMachRate, respectively, and these values are then passed to subsequent functions that 
require these data. The variable MachCasAtTod is a flag that if true, indicates that the Mach-to-CAS 
transitions occurs at the TOD point. 
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If the Mach value for the first waypoint is not set, i.e., the path does not start with a Mach segment, and 
the function terminates with MachCasAtTod set to false. Otherwise, the following is performed. 

if (Crossing Mach first waypoint = 0) terminate this function. Otherwise, 

Set the initial values. 

MachCasAtTod = false 

MachCasModified = false 

Caslndex = index number of the first waypoint 

AltAtMach = 0. 

LastMach = 0 

z = 0 

done = false 

If the TOD Mach data have been modified in a previous invocation of Add Descent Mach Waypoint, 
indicated by a non-empty value for Todld, reset their values. 

if (Todld f empty) 

fini = false 

i = index number of the first waypoint 
Find the waypoint with the name defined in Todld. 
while (fi < (index number of the last waypoint)) and (fini = false )) 
if (Idi = Todld) 
fini = true 

Crossing Mach , = TodMach 
Crossing CASj = 0 
Crossing Rate, = TodMachRate 
Todld = empty string 
i = i + 1 

end of if (Todld f empty) 

Find the first CAS waypoint. 
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fini = false 

i = index number of the first waypoint 

while (fi < index number of the last waypoint) and (fini = false)) 
if (Crossing CAS, > 0) 

Caslndex = i 
fini = true 
i = i + 1 

Determine if the trajectory is already at the CAS altitude, i.e., the initial altitude is the CAS altitude, 
and if so, start in a CAS mode, not Mach. 

// (Crossing Altitudejf rs ( wa yp 0 i n f Altitudec as j nc i ex ) 

done = true 

for (k = index number of the first waypoint; k < Caslndex; k = k + 1) 
if (Crossing Mach k > 0) 

Change the route data so that the trajectory is starting in a CAS mode. 

Invoke the secondary function Mach ToCas. This function is described in a subsequent 
section. 

Crossing CAS k = MachToCas( Crossing Mach k , Altitude caslndex) 

Crossing Mach k = 0 
MachSegment k = false 
end of if (Crossing Mach k > 0) 
if (done = false) 

Find the last Mach value. 
fini = false 

i = index number of the first waypoint 

while ((i < index number of the last waypoint) and (fini = false)) 
if ( Crossing CAS , > 0) fini = true 

else if ( Crossing Mach i > 0) LastMach = Crossing Mach, 
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i = i + 1 


Determine the descent Mach value. 

if (Mach Descent Mach f 0) DescentMach = Mach Descent Mach 

else DescentMach = LastMach 

Determine the Mach-to-CAS transition CAS value. 

if ( Mach Transition CAS > 0) 

MachCas = Mach Transition CAS 

if (Mach Transition CAS < Crossing CAScasindex) MachCas = Crossing CAScasindex 
else MachCas = Crossing CAScasindex 
Find the last Mach altitude. 
fini = false 

i = index number of the first waypoint 

while ((i < index number of the last waypoint) and (fini = false)) 
if (Crossing CAS , > 0) fini = true 

else if (Crossing Altitude, > 0) AltAtMach = Crossing Altitude, 
i = i + 1 

Determine if the Mach is slower than the descent CAS. 

Invoke the secondary function MachCasTransitionAltitude which calculates the altitude where 
the Mach and CAS are equal. This function is described in a subsequent section. 

z = MachCasTransitionAltitude( MachCas, DescentMach ) 

if (z > Crossing Altitudefirst waypoint) 

The path is already below the transition altitude, change the route data so it starts in a CAS 
mode. 

for (k = index number of the first waypoint; k < index number of the last waypoint; k = k + 1) 
done = true 

if (Crossing Machu > 0) 

Crossing CAS) = MachCas 
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Crossing Machk = 0 


MachSegmentk = false 
end of if (done = false) 
if (done = false) 

Find the last Mach value. 
fini = false 

i = index number of the first waypoint 

while ((i < index number of the last waypoint) and (fini = false)) 
if (Crossing CAS , > 0) fini = true 

else if (Crossing Mach, > 0) LastMach = Crossing Mach, 
i = i + 1 

Determine the descent Mach. 

if (Mach Descent Mach f 0) DescentMach = Mach Descent Mach 

else DescentMach = LastMach 

Find the Mach-to-CAS transition CAS. 

if (Mach Transition CAS > 0) MachCas = Mach Transition CAS 
Make sure that the crossing restriction can be obtained. 

if (Mach Transition CAS < Crossing CAS Ca sindex) MachCas = Crossing CAScasindex 
else MachCas = Crossing CAScasindex 
Find the last Mach altitude. 
fini = false 

i = index number of the first waypoint 

while ((i < index number of the last waypoint) and (fini = false)) 
if (Crossing CAS, > 0) fini = true 

else if (Crossing Altitude, > 0) AltAtMach = Crossing Altitude, 
i = i + 1 
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Determine if the Mach is slower than the descent CAS. 


z = MachCasTransitionAltitude( MachCas, DescentMach ) 
if (z > Crossing Altitude firstwaypoint ) 

The path is already below the transition altitude, change the route data so it is starting in a 
CAS mode. 

for (k = index number of the first waypoint; k < index number of the last waypoint; k = k + 1) 
done = true 

if (Crossing Mach k > 0) 

Crossing CAS}- = MachCas 
Crossing Mach k = 0 
MachSegmentk = false 
end of if (done = false) 

If the path still starts with a Mach segment, which may have already been modified in this function, 
test for other special cases. 

if (done = false) 

If required, handle the special case of an accelerated descent. 
if (DescentMach > LastMach) 

Invoke the secondary function HandleDescentAccelDecel. This function handles the special 
case of a Mach acceleration in the descent where the first CAS crossing restriction cannot be 
met. This function is described in a subsequent section. This function may modify the 
waypoint data. 

HandleDescentAccelDecel( Caslndex, LastMach, MachCasModified, DescentMach, 
MachCas ) 

If the descent data are changed, recalculate z. 
if (MachCasModified) 

z = MachCasTransition Altitude ( MachCas, DescentMach ) 

Next, update the waypoint data. 

Mach Descent Mach = DescentMach 
Mach Transition CAS = MachCas 
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end of if (DescentMach > LastMach) 
if (z < Crossing Altitudecasindex) 

At this point, the descent CAS or Mach needs to be changed. 
m = CasToMach( MachCas, Crossing Altitudecasindex) 
if (m > DescentMach) 

Change the descent CAS. 

MachCas = MachToCas( DescentMach, Crossing Altitudecasindex) 

else 

DescentMach = CasToMach( MachCas, Crossing Altitudecasindex) 

Mach Descent Mach = DescentMach 
z = Crossing Altitude casindex 

Perform an extreme limits test, assuming that a valid Mach value will be between 0.6 and 0.9 
Mach. 

if ((DescentMach > 0.9) or (DescentMach < 0.6)) mark this as an error condition 
end of if (z < Crossing Altitude casindex) 

Make sure that there is sufficient distance to slow from the Mach-to-CAS transition speed to 
make the crossing CAS. 

if ( (z> Crossing Altitudecasindex) and (MachCas > Crossing CASc as index) and 
(Crossing Ratecasindex > 0) and (MachCasModified = false)) 

Find the distance atz. This is an iterative solution. 

i = Casindex - 1 

fini = false 

Calculate the headwind at the end point. This calculation the secondary function 
InterpolateWindWptAltitude, described in a subsequent section. 

InterpolateWindWptAltitude( Wind Profilecasindex, Altitudecasindex, Ws, Wd ) 

Headwind = Ws * cosine( Wd - GndTrack casindex) 

CurrentGs = ComputeGndSpeedUsingTrack( Crossing CAScasindex, GndTrackcosindex, 
Altitudecasindex, Ws, Wd) 

Iterate = false 
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OnePass = true 


MCasHold = MachCas 
LastCut = 0 
while (fini = false) 
i = Caslndex - 1 

while ((i > index number of the first waypoint) and (Altitude , < z)) i = i - 1 

if (( Altitude ,• - Altitude i+ i) <0) a = 0 

else a = (z - Altitude i+ i) / (Altitude! - Altitude i+ i) 

Calculate the distance, dx, required to reach the altitude. 
dx = a * ( DTGt - DTG i+ i) + DTG i+ i - DTGcasindex 
InterpolateWindWptAltitude( Wind Profilecasindex, z , Ws2, Wd2 ) 

Hw2 = Ws2 * cosine( Wd2 - GndTrack,) 

AvgHw = (Headwind + Hw2) / 2 

Invoke the secondary function EstimateNextCas. EstimateNextCas is an iterative function 
to estimate the CAS value at the next waypoint. 

CasTest =EstimateNextCas( Crossing CAS Cas index, CurrentGs, true, MCasHold, AvgHw, 
z, dx, Crossing Rate Cas index) 

If required, set up the iteration values, where the iteration value is in CAS. 
if (OnePass = true) 

if (CasTest < MachCas) Iterate = true 

else fini = true 

OnePass = false 

Calculate the iteration step size. 

LastCut = \MachCas - CasTest j 
Limit the step size to no smaller than 2 kt. 
if (LastCut < 2) LastCut = 2 
if (Iterate) 
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if (MachCas > CasTest) s = MachCas - LastCut 
else s = MachCas + LastCut 
LastCut = 0.5 * LastCut 
if (s > MCasHold) s = MCasHold 
Determine if the Mach-to-CAS estimate is valid. 
if (((s + 0.25) > MachCas) and - MachCas\ < 1)) 
fini = true 

Calculate the Mach-to-CAS altitude for the current estimate, 
z = MachCasTransition Altitude ( MachCas, DescentMach ) 

Determine if a deceleration is needed prior to the TOD. Add a 50 ft buffer value. 
if (z > (AltAtMach + 50)) 

Find the TOD waypoint. 
fini2 = false 

j = index number of the first waypoint 

while ((j < index number of the last waypoint) and (ftni2 = false)) 
if (Waypointj is marked as the TOD point) fini2 = true 
else j =j + 1 

The altitude index for the test is the TOD altitude point. 
if (fini2 and (i =j) ) 

Mach Descent Mach = CasToMach( Mach Transition CAS, AltAtMach ) 
MachCas At Tod = true 
end of if (z > (AltAtMach + 50)) 
end of if (((s + 0.25) > MachCas) and (|s - MachCas\ < 1)) 
else 

Mach Transition CAS = s 
MachCas = s 
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z = MachCasTransitionAltitude( MachCas, DescentMach ) 
if (z > Altitude ,) z = Altitudei 

j=j + l 

Add a test to limit the number of iterations to 10. 
if 0 — 10) fini = true 
end of if (Iterate) 
end of while (fini = false) 
end of if (done = false) 

Add Descent Mach Waypoint 

The Add Descent Mach Waypoint function changes the descent waypoint Mach if the descent Mach, 
Mach Descent Mach, is different than the cruise Mach. The function also will add any required, additional 
TCPs. 

The calling program provides as input and retains the subsequent outputs for the following variables: 
Todld, TodMach, and TodMachRate. The variable Todld is the name of the top-of-descent waypoint and 
is initialized as a null string by the calling program. Since this function may overwrite the Mach and 
speed change rate for an input waypoint, these variables allow the function to retain the original values for 
Mach and speed change rate and to then reset these variables to their original values prior to recalculating 
new values. 

If the Mach value for the first waypoint is not set, i.e., the path does not start with a Mach segment, or 
there is no defined descent Mach, i.e., Mach Descent Mach = 0, the function terminates. Otherwise, 

If the previous TOD data for an input waypoint have been changed, these data are restored to their 
original values. 

fini = false 

i = index number of the first waypoint 
The last designated Mach waypoint, 

LastMachlndex = index number of the first waypoint 
The first designated CAS waypoint, 

FirstCasIndex = index number of the first waypoint 
Todlndex = 0 

Find the Mach and CAS waypoints. 
fini = false 
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i = index number of the first waypoint 

while ( (i < index number of the last waypoint) and (fini = false) ) 
if (Crossing Machj > 0) LastMachlndex = i 
else if (Crossing CAS , > 0) 

FirstCasIndex = i 
fini = true 
i = i + 1 

Find the TOD waypoint and Mach. 
fini = false 

i = index number of the first waypoint 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if ((Altitude j < Altitude fi rstwavpoint ) or (Cas Cross , > 0)) 
if (Altitude j f Altitudefirst W ay p0 int) Todlndex = i - 1 
else Todlndex = i 
fini = true 

else if (Crossing Machj > 0) MachAtTod = Crossing Machj 
i = i + 1 

If the vertical segment type has not been defined, mark this as the TOD. 

if ((Todlndex > 0) and (VSegType To didx ~ no type)) VSegTyperodidx = TOD ALTITUDE 

Check for errors. There cannot be a programmed descent Mach if there is a downstream Mach 
restriction. 

if ( (LastMachlndex > Todlndex) or (FirstCasIndex < Todlndex) ) mark this as an error condition 
else 

Save the Mach values for all input waypoints so that they may be reset on subsequent passes back 
to their original input values. 

if (Waypoint Todlndex = input waypoint) 

copy the name of Waypoint Todlndex into Todld 
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TodMach = Crossing Mach To dindex 
TodMachRate = Crossing RateTodindex 
if ( (Waypoint Todindex = input waypoint) and (Crossing Rate To dindex > 0) ) 

CAS Rate = Crossing Rate To dindex 
else CAS Rate = 0. 75 kt / sec (a default value) 

The following is added to force a subsequent speed calculation. 

Crossing Rate To dindex = CAS Rate 

If the aircraft will slow during the descent, do the following: 
if (MachAtTod >Mach Descent Mach) 

Overwrite the TOD Mach value. 

Crossing MachTodindex = Mach Descent Mach 

else 

This is a special case where the aircraft is accelerating to the descent Mach. 

Invoke the secondary function DoTodAcceleration. This function is described in a subsequent 
section. 

DoTodAcceleration( Todldx, MachAtTod ) 

Crossing MachTodindex = MachAtTod 

Compute Mach-to-CAS TCP 

If a Mach-to-CAS transition is required, this functions computes the Mach-to-CAS altitude and inserts 
a Mach-to-CAS TCP. This function is only performed if the input data starts with a Mach Crossing Speed 
for the first waypoint. The function determines the appropriate Mach and CAS values, calculates the 
altitude that these values are equal, and then determines the along-path distance where this altitude occurs 
on the profile. A Mach-to-CAS TCP is then inserted into the TCP list. 

Find the last Crossing Mach and the first Crossing CAS in the list. 

First CAS = 0 

i = index number of the first waypoint 

while ( (i <index number of the last waypoint) and (First CAS = 0) ) 
if (Crossing Mach , > 0) 

Last Mach = Crossing Mach, 
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Last Mach Altitude = Altitude, 
else if (Crossing CAS , > 0) 

First CAS = Crossing CAS , 

CAS Rate = CAS Rate t 
i = i + 1 

If there is a Mach-to-CAS CAS transition speed input, use this value for the First CAS value. 
if (Mach Transition CAS > 0) First CAS = Mach Transition CAS 
Compute the Mach-to-CAS transition altitude, 
z = ComputeMachCasAltitude( FirstCas, LastMach ) 

For an actual implementation, it would be beneficial to check for an error at this point. If z is greater 
than the altitude associated with the Last Mach TCP or if z is less than the altitude associated with the 
First CAS TCP, then an error should be noted. 

Find where z first occurs. 

i = index number of the first waypoint + 1 

finished = false 

while ( (i < index number of the last waypoint) and (finished = false)) 
if (Altitudei > z) i = i + 1 
else finished = true 
Find the distance to this altitude. 
x = Altitude,.! - Altitude , 
if (x < 0) ratio = 0 
else ratio = (z - Altitude,) / x 
d = ratio * (DTG,.i - DTG) + DTG, 

Compute the ground track at distance d along the trajectory and save it as Saved Ground Track. 

Saved Ground Track = GetTrajGndTrk( d ) 

Insert a new TCP at location i in the TCP list. The TCP is inserted between TCP,./ and TCP, from the 
original list. The function InsertWaypoint should be appropriate for the actual data structure 
implementation of this function. 

InsertWaypoint( i ) 
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Mark this TCP as the Mach-to-CAS transition TCP. 


Add the data for this new TCP. 

Crossing Machi = Last Mach 
Crossing CAS , = First CAS 
CAS Ratei = CAS Rate 
DTGj = d 
Altitude , = z 

Crossing Angle, = Crossing Angle i+ i 
Ground Tracki = Saved Ground Track 
Mach , = Last Mach 
CASi = First CAS 

Compute and add the wind data at distance d along the path to the data of TCP,. 
GenerateWptWindProfile( DTG„ TCPi) 

Mark all TCPs from the first TCP (TCPf irs , wa y poim ) to TCP,.] as Mach TCPs. 

Compute Altitude / CAS Restriction TCP 

If an altitude / CAS restriction is required, the Compute Altitude / CAS Restriction TCP function 
computes the altitude / CAS restriction point and insert an altitude / CAS TCP. This is the (U.S.) point 
where the trajectory transitions through 10,000 ft and a 250 kt restriction is required. This function is only 
performed if the previously computed flag NeedlOKRestriction is true. The function determines the 
along-path distance where this altitude / CAS occurs on the profile. A TCP is then inserted into the TCP 
list at this point. The restriction values are Descent Crossing Altitude and Descent Crossing CAS. 

Find the first TCP that is below the Descent Crossing Altitude in the list. 

i = index number of the first waypoint 

k = i 

fini = false 

while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Altitude, < Descent Crossing Altitude) 
k = i 

fini = true 
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i = i + 1 


Find the last CAS restriction prior to the first waypoint below Descent Crossing Altitude, 
i = k - 1 
fini = false 
Last CAS = 0 

while ( (i > 0) and (fini = false) ) 
if (Crossing CAS , > 0) 

Last CAS = Crossing CAS , 
fini = true 
i = i - 1 

Determine if an altitude or CAS TCP is required. If it is, add it. 
if ( (TCP k is a Mach segment) and (Last CAS > Descent Crossing CAS) ) 
i = k; 

Find the distance to this altitude. 
x = Altitude,.! - Altitude ,• 
if (x < 0) ratio = 0 

else ratio = (Descent Crossing Altitude - Altitude ) / x 
d = ratio * (DTGj_! - DTGj + DTGj 

Compute the ground track at distance d along the trajectory and save it as Saved Ground Track. 
Saved Ground Track = GetTrajGndTrk( d ) 

Insert a new TCP at location i in the TCP list. The TCP is inserted between TCP,./ and TCP, from 
the original list. The function InsertWaypoint should be appropriate for the actual data structure 
implementation of this function. 

InsertWaypointf i ) 

Mark this TCP as the altitude / CAS restriction TCP. 

VSegType , = altitude CAS restriction 
Turn Type, = no turn 
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Add the data for this new TCP. 

Crossing Mach , = 0 

Crossing CAS t = Descent Crossing CAS 
Use a high value, arbitrary CAS rate. 

CAS Ratei = 0. 75 kt / sec 
DTG, = d 

Altitude i = Descent Crossing Altitude 
Crossing Angle, = Crossing Angle i+ i 
Set the Mach flag for TCP, to false 
Ground Trackj = Saved Ground Track 
Mach i = 0 

CASj = Descent Crossing CAS 

Compute and add the wind data at distance d along the path to the data of TCP ,. 
GenerateWptWindProfile( DTG „ TCP,) 

Test for Altitude / CAS Restriction Requirement 

The Test for Altitude / CAS Restriction Requirement function determines if the addition of an altitude / 
CAS restriction point is required. This is the (U.S.) point where the trajectory transitions through 10,000 
ft and a 250 kt restriction is required. This function determines the value of the NeedlOKRestriction flag. 
The function can only be called after an initial, preliminary trajectory has been generated. The restriction 
values are Descent Crossing Altitude and Descent Crossing CAS. 

NeedlOKRestriction = false 

if ( (Descent Crossing Altitude > 0) and (Descent Crossing CAS > 0) ) ok = true 
else ok = false 

If we don't start above 10,000ft, skip this whole routine. 
if (ok and ( Altitude f irstwaypoint > Descent Crossing Altitude)) 

Find the first point below Descent Crossing Altitude 
fini = false 
i = 0 
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while ( (i <index number of the last waypoint) and (fini = false) ) 
if (Altitudei < Descent Crossing Altitude) 

Find the distance to this altitude. 
x = Altitude,.! - Altitude ,■ 
if (x < 0) ratio = 0 

else ratio = (Descent Crossing Altitude - Altitudei) / x 
s = ratio * (CAStr! - CAS) + CAS, 

if (s > (Descent Crossing Cas + 2) ) NeedlOKRestriction = true 
fini = true 
i = i + 1 

Compute TCP Speeds 

The Compute TCP Speeds function is similar to Compute TCP Altitudes in its design. Beginning with 
the last waypoint, this function computes the Mach or CAS at each previous TCP and inserts any 
additional speed TCPs that may be required to denote a change in the speed profile. The function uses the 
current speed constraint, searches backward for the previous constraint, and then computes the distance 
required to meet this previous constraint. The speeds for all of the TCPs within this distance are computed 
and added to the data for the TCPs. If the along-path distance to meet the previous constraint is not at a 
TCP, a new speed VTCP is inserted at this distance. This function invokes two secondary functions, 
described in the subsequent text, with the invocation dependent on the constraint speed, whether it is a 
Mach or a CAS value. This function is performed in the following steps: 

Set the current constraint index number, cc, equal to the index number of the last waypoint, 
cc = index number of the last waypoint 

The speed of the first waypoint is set to its crossing speed. 

if (Crossing Mach firstwaypoint > 0) 

Much Ji rs t waypoint = Crossing Mach firs , way point 

CASfirst waypoint MachToC as( Mach j, rs , waypoint. Altitude f irs t waypoint) 


else 


CAS f, rs t waypoint Crossing T AS first waypoint 


Mach frst waypoint CasToMach( C AS first waypoint, Altitude first waypoint) 

The speed of the last waypoint is set to its crossing speed. 


CAS cc = Crossing CAS cc . 
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A flag signifying that Mach segment computation has begun is set to false, 

Doing Mach = false 

While (cc > index number of the first waypoint) 

Set the Mach flag if the current TCP is the Mach-to-CAS transition point. 
if (TCP cc = Mach Transition CAS) Doing Mach = true 
if (Doing Mach) ComputeTcpMach(cc) 
else ComputeTcpCas(cc) 

end of while cc > index number of the first waypoint 

Compute Secondary Speeds 

The Compute Secondary Speeds function adds the Mach values to CAS TCPs, the CAS values to Mach 
TCPs, and the ground speed values to all TCPs. This function is performed in the following steps: 

Doing Mach = false 

Working backwards from the runway, compute the relevant speeds. 

for (i = index number of the last waypoint; i > index number of the first waypoint; i = i - 1) 

Set the flag if the current TCP is the Mach-to-CAS transition point. 
if (TCP, = Mach Transition CAS) Doing Mach = true 
if (Doing Mach) Cas, = MachToCas( Mach it Altitudei) 
else Mach i = CasToMach( Cas i} Altitudei) 

Compute the ground track. 

if (i = index number of the first waypoint) track = Ground Track , 
else if (WptlnTurn(i) or (TCP, = turn-exit)) track = Ground Tracks 
else track = Ground Tracks/ 

Compute the ground speed. This also requires the computation of the wind at this point. 

InterpolateWindWptAltitude( Wind Profile b Altitudei, Wind Speed, Wind Direction ) 

Ground Speed , = ComputeGndSpeedUsingTrack ( Cos,, track, Altitudei, Wind Speed, 

Wind Direction ) 

end of for (i = index number of the last waypoint; i > index number of the first waypoint; i = i - 1) 
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Update Turn Data 

The Update Turn Data function computes the turn data for each turn waypoint and modifies the 
associated waypoint's turn data sub-record. This function performs as follows: 

KtsToFps = 1.69 

Nominal Bank Angle = 22 

index = index number of the first waypoint + 1 

while (index < index number of the last waypoint) 

Find the next input waypoint with a turn. 

while ( (index < index number of the last waypoint) and ( (TCP index f input waypoint) or 
(not WptInTurn( index )))) index = index + 1 

If there are no errors and there is a turn of more than 3 -degrees, compute the turn data. 

if (index < index number of the last waypoint) 

Find the start of the turn. 

i = index - 1 

while (TCP f turn-entry) i = i - 1 
start = i 

The following are all approximations and are based on a general, constant radius turn. 

The start of turn to the midpoint data is as follows, noting that the ground speeds for all points 
must be valid at this point. 

The overall distance d for this part of the turn is, 

d = DTG start - DTGindex 

The special case with 0 distance between the points is, 

if (d < 0) AvgGsFirstHalf = (Ground Speeds, art + Ground Speed im u, x ) / 2 

else 

The overall average ground speed is computed as follows, noting that it is the sum of 
segment distance / overall distance * average segment ground speed. 

AvgGsFirstHalf = 0 

for (j = start; j < (index - 1); j =j + 1) 

dx = DTGj - DTGj+i 
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AvgGsFirstHalf = AvgGsFirstHalf + (dx / d) 

* (Ground Speed ) + Ground Speedy f) / 2 

Now, find the end of the turn. 

i = index + 1 

while (TCPi f turn-exit) i = i + 1 
end = i 

Now, find the midpoint to the end of the turn. 

The overall distance for this part of the turn is, 

d = D TG index - DTGend 

Test for the special case, 0 distance between the points. 
if(d<0) 

AvgGsLastHalf = (Ground Speed ,„* x + Ground Speed er f / 2 

else 

Compute the overall average ground speed noting that it is the sum of segment 
distance / overall distance * average segment ground speed. 

AvgGsLastHalf = 0 

for (j = index; j < (end - 1); j =j+ 1) 

dx =DTGj - DTGj+i 

AvgGsLastHalf = AvgGsLastHalf + (dx / d) * 

(Ground Speed j + Ground Speed )+i) / 2 

end of for (j = index; j < (end - 1); j =j + 1) 

end of else if (d < 0) 

full turn = DeltaAngle(Ground Track starb Ground Track en f 
half turn =full turn / 2 

Compute the outputs from the average ground speed. 

Average Ground Speed = (AvgGsFirstHalf + AvgGsLastHalf / 2 
Save the ground speed data in the turn data for this waypoint. 
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Turn Data Average Ground Speed index = Average Ground Speed 

Compute the turn radius and associated data. This set of calculations is not performed if the 
waypoint is a special, RF center-of-tum turn waypoint. 

if (Wpti f rf-turn-center) 

The general equation is turn rate = c tan(bank angle) / v. If the bank angle is a constant, 
turn rate = cO / v. The Nominal Bank Angle = 22 degrees. 

cO = 57.3 * 32.2 / KtsToFps * tangent(Nominal Bank Angle) 

w = cO / Average Ground Speed 

The time to make the turn is, 

Turn Data Turn Time index = \full turn \ / w 

The turn radius is, 

Turn Data Turn Radius in de X = (57.3 * KtsToFps * Average Ground Speed) / 

(6076 *w) 

The along-path distance for the turn is, 

Turn Data Path Distance in d ex = {full turn\ * Turn Data Turn Radius 57.3 

else 

These are the data for an RF turn. The along-path distance for the turn is, 

Turn Data Path Distance in d ex = full turn\ * Turn Data Turn Radius in d ex / 57.3 
The time to make the turn is, 

Turn Data Turn Time in d ex = Turn Data Path Distance in d ex / Average Ground Speed 

* 3600 

Save the turn data for the first half of the turn, denoted by the " 1 " in the variable name. 

Turn Data Casl index = CAS slarl 

Turn Data Average Ground Speedl m d ex = AvgGsFirstHalf 
Turn Data Trackl in d ex = Ground Track star t 

The Straight Distance values are the distances from the turn-entry TCP to the waypoint and 
from the waypoint to the turn-exit TCP. See the example in figure 6. 

Turn Data Straight Distance 1 in d ex = Turn Data Turn Radius in d ex * tangent(\half turn \) 
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Figure 6. Turn distances for waypoint,. 

The Path Distance values are the along-the-path distances from the turn-entry TCP to a point 
one-half way along the turn and from this point to the turn-exit TCP. See the example in 
figure 6. 

Turn Data Path Distancel in d ex = | half turn\ * Turn Data Turn Radius index / 57.3 

Compute the midpoint waypoint data. This set of calculations is not performed if the 
waypoint is a special, RF center-of-tum waypoint. 

if (Wpti f rf-turn-center) 

w = cO / AvgGsFirstHalf 

Turn Data Turn Timel index = \half turf / w 


else 


These are the data for an RF turn. 

Turn Data Turn Timel in dex — Turn Data Path Distancel index / AvgGsFirstHalf * 3600 

The data for the midpoint to the end of the turn, denoted by the "2" in the variable name, are 
as follows: 

Turn Data Cas2 ind e X = CAS en d 

Turn Data Average Ground Speed2 index = AvgGsLastHalf 
Turn Data Track2 index = Ground Track end 

The distances for the second half of the turn are the same as for the first, but their calculations 
are recomputed here for clarity. 

Turn Data Straight Distance2 index = Turn Data Turn Radius in d ex * tangent(\half turf) 

Turn Data Path Distance2 index = \half turf * Turn Data Turn Radius in d ex / 57.3 

Compute the data for the last half of the turn. Again, this set of calculations is not performed 
if the waypoint is a special, RF center-of-tum waypoint. 
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if (Wpti f rf-turn-center) 
w = cO / AvgGsLastHalf 
Turn Data Turn Time2 index = | half turn\ / w 


else 


These are the data for an RF turn. 

Turn Data Turn Time2 index = Turn Data Path Distance2 index / AvgGsLastHalf * 3600 
The DTG values are as follows: 

DTG start = DTG i ndex + Turn Data Path Distance 1 index 
DTG end = DTG index - Turn Data Path Distance2 index 

Since the turn waypoints have been moved, the wind data need to be updated for the new 
locations. 

if (TCP start T input waypoint) GenerateWptWindProfile( DTG sta rt, TCP sta rt) 
if (TCP end f input waypoint) GenerateWptWindProfile( DTG end , TCP end ) 
end of if (index < index number of the last waypoint) 
index = index + 1 

end of while (index < index number of the last waypoint) 

Delete TCPs 

The Delete TCPs function deletes the altitude, speed, and Mach-to-CAS TCPs. The remaining TCPs 
will only consist of input waypoints, turn-entry, and turn-exit TCPS. This function also removes any flags 
that associate any remaining TCPs with a speed or altitude change, e.g., a waypoint marked as the 10,000 
ft, 250 kt restriction. 

Update DTG Data 

The Update DTG Data function is performed after the turn data have been updated and the VTCPs 
have been deleted. Only input, turn-entry, and turn-exit TCPs should be in the list at this time. If the input 
test flag, TestOnly, is true, then only the testing portions of this function are used. 

if (TestOnly = false) DTG first way point = 0 

i = index number of the last waypoint 

while (i > index number of the first waypoint) 

Determine if there is a turn at either end and adjust accordingly. 

if (WptInTurn( i )) 
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if (TestOnly = false) DTG,_i = DTG , + Turn Data Path Distance 1 , 

The following is the difference between going directly from the waypoint to going along the 
curved path. 

PriorDistanceOffset = Turn Data Straight Distance 1 , - Turn Data Path Distance I , 
else PriorDistanceOffset = 0 
Find the next input waypoint. 
n = i - 1 

while (TCP n f input waypoint) n = n - 1 
if (WptInTurn( n )) 

The following is the difference between going directly from the waypoint to going along the 
curved path. 

DistanceOffset = Turn Data Straight Distance2„ - TurnData.PathDistance2„ 

The DTG to the input waypoint is then: 

if (TestOnly = false) DTG n = (Center to Center Distance , - PriorDistanceOffset - 
DistanceOffset) + DTG, 

If the DistanceOffset is greater than Center to Center Distance,, then the turn is too big. 
if (DistanceOffset > Center to Center Distance ) mark this as an error condition 
The turn-exit DTG is then, 

if (TestOnly = false) DTG n +i = DTG n - Turn Data Path Distance2„ 
else if (TestOnly = false) 

The next waypoint is not in a turn. 

DTG n = Center to Center Distance , - PriorDistanceOffset + DTG , 
i = n 

end of while (i > 0) 
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Check Turn Validity 

The Check Turn Validity function is performed after the turn data have been updated and the VTCPs 
have been deleted. Only input, turn-entry, and turn-exit TCPs should be in the list at this time. The 
function simple checks that there are no turns within turns by examining the DTG values. 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 

if (DTG, < DTG i+I ) mark this as an error condition 

Recover the Initial Mach Segments 

This function. Recover the Initial Mach Segments, attempts to recover the Mach portion of the 
trajectory if the initial segments should be Mach but have been internally converted to CAS in the 
function Meet Cruise CAS Waypoint Restriction. This function uses the Mach value that was saved at the 
start of this program from the first waypoint of the original route. This saved Mach value, First Waypoint 
Mach, is compared to the Mach equivalent value of the CAS at the initial waypoints and if these Mach 
values are the same, these waypoints are marked as Mach segments instead of CAS segments. 

Only perform this function if the calculated trajectory does not start with a Mach segment but the 
original route does start with a Mach value. 


if ((Mach Segment index number of the firs, waypoint = false) and (First Waypoint Mach fO)) 


Mach CasToMach( Cl ossing C AS index number of the first waypoint* 


Altitude index number of the first waypoint) 


if (Mach ~ First Waypoint Mach) 


fini = false 

i = index number of the last waypoint 


FirstCas = Crossing CAS, 


index number of the first waypoint 


If there is no Mach transition altitude set, set the transition values. 


if (Mach Transition Altitude = 0) 

Mach Descent Mach = First Waypoint Mach 
Mach Transition Mach = First Waypoint Mach 
Mach Transition Cas = FirstCas 


Mach Transition Altitude = Altitude in dex of first waypoint 

while (( i < (index number of the last waypoint - 1)) and (fini = false)) 

Test that the CAS computed for the waypoint is the same as the FirstCas, that except for 
the first waypoint that there is not speed crossing condition at the waypoint, and that the 
altitude computed for the waypoint is the same as the altitude for the first waypoint. 

if (( Cast = FirstCas) and ( (i = index number of the last waypoint) or 
((Crossing Mach i = 0) and ( Crossing CASi = 0))) and 
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(Altitude j Cl OSS il l g Altitude index number of the first waypoint)) 

If the previous conditions are turn, set this waypoint as a Mach segment. 

Mach Segment, = true 

Change the speed crossing values for the first waypoint. 
if ( Crossing CAS , > 0) 

Crossing CAS , = 0 

Crossing Mach, = First Waypoint Mach 
end of if (( CaSj = FirstCas)...) 
else fini = true 
i = i + 1 

Insert CAS Descent VTCPs 

This function inserts vertical TCPs between constant CAS descent waypoints to improve the TAS 
estimation when using the data provided by this algorithm. This updating occurs at 3,000 ft intervals. 

Update Altitude = 3000 

Find the first CAS point. 

j = 0 

while ( (Mach Segment , = true) and (j < index number of the last waypoint)) j =j + 1 
for ( i = j; i < (index number of the last waypoint - 1); i = i + 1) 

DeltaZ = Altitude, - Altitude, , / 

Update at 3000 ft intervals but skip the update if the waypoint is within 500 ft of the test altitude. 
if ( (DeltaZ > (Update Altitude + 500)) and (Cas, ~ Cas, + /)) 
z = Altitudei - Update Altitude 
dx = DTG, - DTG, + / 
a = arctangent2 (DeltaZ, 6076 * dx) 
d = DTG, - Update Altitude / tan(a) / 6076 

Compute the ground track at distance d along the trajectory and save it as Saved Ground 
Track. 
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Saved Ground Track = GetTrajGndTrk( d ) 
k = i + 1 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted between TCP/,./ and 
TCP, ’ from the original list. The function InsertWaypoint should be appropriate for the actual 
data structure implementation of this function. 

InsertWaypoint( k ) 

Update the waypoint-type data in the new waypoint. 

WptType k = VTCP 
VSegType k = TAS adjustment 
TurnTypek = no turn 

Update the crossing data in the new waypoint. 

Crossing Machk = 0 
Crossing CASk = 0 
Crossing Rate k = 0 
CASk = CAS k+ i 
DTGk = d 
Altitudek = z 

Machk = CasToMach( CASk, Altitudek) 

Mach Segmentk = false 
Crossing Angle k = Crossing Angle k +i 
Ground Trackk = Saved Ground Track 
Compute and add the wind data at this waypoint. 

GenerateWptWindProfile( DTGk, TCPk) 

Compute the wind at the waypoint altitude and then waypoint's ground speed. 

In tei po late Win d Wp tA ltitude( Wind Profde k , Altitudek, Ws, Wd ) 

Ground Speedk = ComputeGndSpeedUsingTrack( CAS k , Ground Trackk-i, Altitudek, Ws, Wd ) 
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Compute TCP Times 

The function Compute TCP Times calculates the time to each TCP. The calculations begin at the 
runway (the last waypoint), working backwards, and compute the TTG to each TCP. 

TTG index number of the last waypoint 0 

for (i = index number of the last waypoint; i > index number of the first waypoint; i = i - 1) 

Average Ground Speed = (Ground Speedy Ground Speed,) / 2 
x = DTGj-i - DTGi 

Test for an error condition where the distance is less than 0. 

fi (x < 0) 

If the distance is close to 0, e.g., within 200 ft, set the distance to the previous and ignore the 
error. 

if (x> (-200 / 6076)) 

DTGi = DTGi., 
x = 0 

Allow a larger margin of error for an RF turn. 

else if ( (x >-0.05) and (TurnType, = turn-entry) ) and (Center Of Turn Latitudeif 0 )) 

DTGi = DTGi., 
x = 0 

else mark this as an error condition 
Delta Time = 3600 * x / Average Ground Speed 
TTG,., = TTG, + Delta Time 

Compute TCP Latitude and Longitude Data 

With the exception of the input waypoints, the Compute TCP Latitude and Longitude Data function 
computes the latitude and longitude data for all of the TCPs. 

In Turn = false 

Last Base = index number of the first waypoint 
Next Input = index number of the first waypoint 
Turn Index = index number of the first waypoint 
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Turn is Clockwise = true 


Turn Adjustment = 0 

Base Latitude = Latitude Last Base 

Base Longitude = Longitude Last Ba se 

for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 
if (TCP, = turn-entry) 

Turn Adjustment = 0 
InTurn = True; 

Find the major waypoint for this turn. 

Next Input = i + 1 

while ( (TCP Next input f input waypoint) and (Next Input < index number of the last waypoint) ) 
Next Input = Next Input + 1 

Turn Index = Next Input 

a = DeltaAnglel Ground Tracks Ground Track^ ex t input) 
x = Turn Data Turn Radius mm index/ cosine(a) 
if (a > 0) Turn Clockwise =true 
else Turn Clockwise = false 

if (Turn Clockwise) al = Ground Trackpum index + 90 
else al = Ground Trackpum index - 90 

Now compute the relative latitude and longitude values. The function RelativeLatLon is 
described in a subsequent section. 

RelativeLatLong(Latitude TurnIndex , Longitude Turn Index , al, x), returning Center Latitude and 
Center Longitude 

end of if (TCPi = turn-entry) 

if (In Turn) 

Turn Adjustment = 0 

if (Turn Clockwise) al = Ground Track , - 90 
else al = Ground Track , + 90 
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if (TCPj = input waypoint) 

Turn Data Center Latitudei = Center Latitude 
Turn Data Center Longitudei = Center Longitude 

RelativeLatLong(C enter Latitude, Center Longitude, al, Turn Data Turn Radius Tuni index), 
returning Turn Data Latitudei and Turn Data Longitude , 

end of if (TCPt = input waypoint) 

else RelativeLatLon(Center Latitude, Center Longitude, al, Turn Data Turn Radius Nex , input), 
returning Latitudei and Longitudei 

if (TCP, = turn-exit) 

Turn Adjustment = Turn Data Straight Distanced TumIndex - 
Turn Data Path Distanced Turn lndex 

In Turn = false 
Last Base = Next Input 
Base Latitude = Latitude Last BaS e 
Base Longitude = Longitude L asl Base 
end of if (In Turn) 
else 

if (TCPj = input waypoint) 

Turn Adjustment = 0 
Last Base = i 

Base Latitude = Latitude i as t Base 
Base Longitude = Longitude Last Base 

else 

RelativeLatLong(Base Latitude, Base Longitude, Ground Tracks, DTGiastnase - DTG, + 
Turn Adjustment), returning Latitudei and Longitudei 

end of for (i = index number of the first waypoint; i < index number of the last waypoint; i = i + 1) 
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Description of Secondary Functions 

The secondary functions are listed in alphabetical order. Note that standard aeronautical functions, such 
as CAS to Mach conversions, CasToMach, are not expanded in this document but may be found 
numerous references, e.g., reference 22. It may also be of interest to include atmospheric temperature or 
temperature deviation in the wind data input and calculate the temperature at the TCP crossing altitudes to 
improve the calculation of the various speed terms. 

BodDecelerationDistance 

The function BodDecelerationDistance estimates the distance required for the special case of a 
deceleration to a CAS restricted waypoint from the Mach-to-CAS transition. This function is invoked 
from HandleDescentAccelDecel, which passes in the index number for the bottom-of-descent (TOD) 
waypoint, Bodlndex, the Mach transition to CAS altitude, MachTransitionAlt, and the CAS at the Mach 
transition to CAS, TransitionCas. The function returns the distance from the index point of the 
deceleration, Distance. 

Estimate the distance to the new Mach value. Begin by finding the time to do the deceleration. 
t = (TransitionCas - Crossing CAS Bo didx) / Crossing Rate Bo didx 
Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude( Wind Profile Bodidx, Altitude Bo didx, Ws, Wd ) 

Calculate the ground track at the current point. 
if (WptInTurn( Bodidx )) track = Ground Track Bo didx-i 
else track = Ground Track Bo didx 
Calculate the ground speed over this segment. 

BodGs = ComputeGndSpeedUsingTrack( Crossing CAS Bo didx, track, Altitude, Bodidx, Ws, Wd ) 

DescentGs = ComputeGndSpeedUsingTrack( TransitionCas, track, MachTransitionAlt, Ws, Wd ) 
Calculate the average groundspeed, AvgGS. 

AvgGs = (BodGs + DescentGs) / 2 
The distance estimate is AvgGs * t . 

Distance = AvgGs * t / 3600 
ComputeGndSpeedUsingMachAndTrack 

The ComputeGndSpeedUsingMachAndTrack function computes a ground speed from track angle 
(versus heading), track, Mach, Mach, altitude, Altitude, and wind data, Wind Speed and Wind Direction. 

CAS = MachToCas( Mach, Altitude ) 

Ground Speed = ComputeGndSpeedUsingTrack( CAS, track, Altitude, Wind Speed, Wind Direction ) 
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ComputeGndSpeedUsingTrack 

The ComputeGndSpeedUsingTrack function computes a ground speed from track angle (versus 
heading), track, CAS, CAS, altitude. Altitude, and wind data, Wind Speed and Wind Direction. 

b = DeltaAngle( track, Wind Direction ) 

if (CAS <0) r = 0 

else r = ( Wind Speed / CasToTas Conversion! CAS, Altitude )) * sine(b) 

Limit the correction to something reasonable. 
if (r\ > 0.8) r = 0.8 * r/\r\ 
heading = track + arcsine(r) 
a = DeltaAngle( heading, Wind Direction ) 

TAS = CasToTas Conversion( CAS, Altitude ) 

Ground Speed = (Wind Speed 2 + TAS 2 - 2 * Wind Speed * TAS * cosine(a) ) 0 ' 5 

ComputeGndTrk 

The ComputeGndTrk function computes the ground track at the along-path distance equal to distance., 
where distance must lie between TCP-,.i and TCP i+ i. It is assumed that the value for Ground Track, is 
invalid. The function uses a linear interpolation based on DTG,./ and DTG i+ i, with the index value i input 
into the function and where the distance, distance, must lie between these points. 

d = DTG i-i - DTG ,+y 

if (d < 0) Ground Track = Ground Tracks 
else 


a = (1 - (distance - DT i+1 ) / d) * DeltaAngle(Ground Track,./, Ground Track i+ i) 

Ground Track = Ground Track,, i + a 

ComputeTcpCas 

The index variable cc is passed into and out of the ComputeTcpCas function. Beginning with the last 
waypoint, this function computes the CAS at each previous TCP and inserts any additional speed TCPs 
that may be required to denote a change in the speed profile. The function uses the current speed 
constraint, searches backward for the previous constraint, and then computes the distance required to meet 
this previous constraint. The speeds for all of the TCPs within this distance are computed and added to the 
data for the TCPs. If the along-path distance to meet the previous constraint is not at a TCP, a new speed 
VTCP is inserted at this distance. Because there is no general closed form solution to compute distances 
to meet the deceleration constraints, an iterative technique is used in this function. This function is 
performed in the following steps: 

While ((cc > index number of the first waypoint) and (TCP cc f Mach Transition CAS)) 

Determine if the previous constraint cannot be met. 
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If (CAS cc > Crossing CAS cc ) 

If this is the last pass through the algorithm, mark this as an error condition 
CAS cc = Crossing CAS cc 

Find the prior waypoint index number pc that has a CAS constraint, e.g., a crossing CAS 
( Crossing CAS pc f 0). This may not always be the previous (i.e., cc - 1) waypoint. 

The initial condition is the previous TCP. 

pc = cc - 1 

while ( (pc > index number of the first waypoint) and (TCP pc f Mach Transition CAS) 
and (Crossing CAS pc = 0) ) pc = pc - 1 

Save the previous crossing speed. 

Prior Speed = Crossing CAS pc 

Save the current crossing speed (Test Speed) at TCP,,,, and the deceleration rate ( Test Rate ) noting 
that the first and last waypoints always have speed constraints and except for the first waypoint, 
all constrained speed points must have deceleration rates. 

Test Speed = Crossing CAS cc 

Test Rate = Crossing Rate cc 

Compute all of the TCP speeds between the current TCP and the previous crossing waypoint. 
k = cc 
while k> pc 

If the previous speed has already been reached, set the remaining TCP speeds to the previous 
speed. 

if (Prior Speed < Test Speed) 

for (k = k- 1; k> pc; k = k - 1) 

CASk = Test Speed 
Machk = CasToMach( CASk, Altitude 
Set the speeds at the last test point. 

CAS P c = Test Speed 

if (Mach P c = 0) Mach pc = CasToMach( CAS pc , Altitude pc ) 

else 
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Estimate the distance required to meet the crossing restriction using the winds at the 
current altitude. This is a first-estimation. 

Compute the time to do the deceleration. 

t = ( Prior Speed - Test Speed) / Test Rate 

Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude( Wind Profile^, Altitude/, Wind Speedl, Wind Directionl ) 
The ground track at the current point is, 
if (WptlnTunf k )) Track = Ground Tracks 
else Track = Ground Tracks/ 

Current Ground Speed = ComputeGndSpeedUsingTrack( Test Speed, Track, 

Altitude k , Wind Speedl, Wind Directionl ) 

Compute the wind speed and direction at the prior altitude. 

InterpolateWindWptAltitude( Wind Profile k _i, Altitude h Wind Speedl, Wind Directionl ) 
The ground speed at the prior point. 

Prior Ground Speed = ComputeGndSpeedUsingTrack( Prior Speed, GndTrack k ./, 
Altitude k .i, Wind Speedl, Wind Directionl ) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2 

The distance estimate, dx, is Average Ground Speed * t. 

dx = Average Ground Speed * t / 3600 

Recalculate the distance required to meet the speed using the previous estimate distance 
dx. 

Begin by computing the altitude, AltD, at distance dx. 

if (Altitude k > Altitude k .i) AltD = Altitude k 

else 

AltD = (6076 * dx) * tangent! Crossing Angle/J + Altitude k 
if (AltD >Altitude k _i) AltD = Altitude k 
The new distance x is DTG k + dx. 

Compute the winds at AltD and distance x. 
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InterpolateWindAtDistance( AltD, x, Wind Speed2, Wind Direction2 ) 

The track angle at this point, with GetTrajGndTrk defined in this section: 

Track2 = GetTrajGndTrk( x ) 

The ground speed at altitude AltD is then, 

Prior Ground Speed = ComputeGndSpeedUsingTrack( Prior Speed, Track2, AltD, 

Wind Speed2, Wind Direction2 ) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2 

dx = Average Ground Speed * t / 3600 

If there is a TCP prior to dx, compute and insert its speed. 

If the distance is very close to the waypoint, just set the speed. 
if ( (DTGk-i < (DTGk + dx + some small value) ) 

if (\ DTGk-i - DTGk - dx \ < some small value) CASk-i = Prior Speed 
else 

Compute the speed at the waypoint using v 2 = Vo 2 + 2ax to get v. 

The headwind at the end point is, 

HeadWind2 = Wind Speed2 * cosine(Wind Direction2 - Ground Track k .i) 
dx = DTGk-i - DTGk 

The value of CASk-i is computed using function EstimateNextCas, described in 
this section. 

CASk-i = EstimateNextCas ( Test Speed, Current Ground Speed, false, 

Prior Speed, Head Wind2, Altitudek, dx, Crossing Rate cc ) 

Determine if the constraint is met. 

tf ( (k-1) = pc) 

Determine the allowable crossing window, accounting for special conditions. 

if (( (pc + 1) < index number of the last waypoint) and 
(VSegType pc = MACH CAS) ) CrossingWindow = 5 

else CrossingWindow = 1 

Was the crossing window speed met? If not, set this as an error. 
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if (\ CAS pc - Crossing CAS pc \ > CrossingWindow) 

Mark this as an error condition 

Always set the crossing exactly to the crossing speed. 

CAS pc = Crossing CAS pc 

Set the test speed to the computed speed. 

Test Speed = CAS k .i 

Back up the index counter to the next intermediate TCP. 
k = k- 1 

end of if ( (DTG k .i < (DTG k + dx + some small value) ) 
else 

The constraint occurs between this TCP and the previous TCP. A new VTCP needs 
to be added at this point. 

The along path distance d where the VTCP is to be inserted is: 
d = DTG k + dx 

Save the ground track value at this distance. 

Saved Ground Track = GetTrajGndTrk( d ) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted between 
TCP/,./ and TCP* from the original list. The function InsertWaypoint should be 
appropriate for the actual data structure implementation of this function. 

InsertWaypoint( k ) 

Update the data for the new VTCP which is now TCP k . 

WptType k = VTCP 

if (VSegType k = no type) VSegType k = SPEED 
TurnType k = no turn 
DTG k ~ d 

The altitude at this point is computed as follows, recalling that the new waypoint is 
TCP k : 

if (Altitude k +i >Altitude k -i) Altitude k = Altitude k -i 
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else Altitude^ = (6076 * dx) * tangent(Crossing Angle^+i) + Altitude k+i 
CASk = Prior Speed 

Add the ground track data which must be computed if the new VTCP occurs within a 
turn. The functions WptlnTurn and ComputeGndTrk are described in this sections. 

if (WptInTurn( k )) Ground Tracks = ComputeGndTrk( k, d ) 

else Ground Tracks = Saved Ground Track 

Compute and add the wind data at distance d along the path to the data of TCPk ■ 
GenerateWptWindProfile( d, TCPk) 

Test Speed = Prior Speed 

Since TCPk, has now been added prior to pc, the current constraint counter cc needs 
to be incremented by 1 to maintain its correct position in the list. 

cc = cc + 1 

end of while k > pc. 

Now go to the next altitude change segment on the profile. 
cc = k 

end of while cc > index number of the first waypoint 

ComputeTcpMach 

The index variable cc is passed into and out of the ComputeTcpMach function. This function is similar 
to ComputeTcpCas with the exception that the computed Mach rate will need to be recomputed with any 
change of altitude. Beginning with the last Mach waypoint (the Mach waypoint that is closest to the 
runway in terms of DTG), this function computes the Mach at each previous TCP and inserts any 
additional speed TCPs that may be required to denote a change in the speed profile. The function uses the 
current speed constraint, searches backward for the previous constraint, and then computes the distance 
required to meet this previous constraint. The speeds for all of the TCPs within this distance are computed 
and added to the data for the TCPs. If the along-path distance to meet the previous constraint is not at a 
TCP, a new speed VTCP is inserted at this distance. Because there is no general closed form solution to 
compute distances to meet the deceleration constraints, an iterative technique is used in this function. This 
function is performed in the following steps: 

While (cc > index number of the first waypoint) 

Determine if the previous constraint cannot be met. 

If (Mach cc > Crossing Mach cc ) 

If this is the last pass through the algorithm, mark this as an error condition 
Machec = Crossing Mach cc 


55 



Find the prior waypoint index number pc that has a Mach constraint, e.g., a crossing Mach 
(■ Crossing Mach pc f 0). This may not always be the previous (i.e., cc - 1 ) waypoint. 

Initial condition is the previous TCP. 

pc = cc - 1 
finished = false 

while ( (pc > index number of the first waypoint) and (TCP pc f Mach Transition CAS) 
and (Crossing CAS pc = 0) ) pc = pc - 1 

Save the previous crossing speed. 

Prior Speed = Crossing Mach pc 

Save the current crossing speed {Test Speed) at TCP,, and the deceleration rate ( Test Rate) noting 
that the first and last waypoints always have speed constraints and except for the first waypoint, 
all constrained speed points must have deceleration rates. 

Test Speed = Crossing Mach cc 

Test Rate = CasToMach( Altitude cc , Crossing Rate cc ) 

Compute all of the TCP speeds between the current TCP and the previous crossing waypoint. 
k = cc 
while k> pc 

If the previous speed has already been reached, set the remaining TCP speeds to the previous 
speed. 

if (Prior Speed < Test Speed) 

for (k = k- 1; k> pc; k = k - 1) 

Machk = Test Speed 

CAS k = MachToCas( Mach k , Altitude k ) 

Mark TCP k as a Mach segment. 

Set the speeds at the last test point. 

Mach pc = Test Speed 

CAS pc = MachToCas( Mach pc , Altitude pc ) 
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Estimate the distance required to meet the crossing restriction using the winds at the 
current altitude. This is a first-estimation. 

Compute the time to do the deceleration. 

t = ( Prior Speed - Test Speed) / Test Rate 

Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude( Wind Profile^ Altitudek, Wind Speedl, Wind Directionl ) 
The ground track at the current point is, 
if (WptlnTurn(k)) Track = Ground Tracks 
else Track = Ground Tracks i 

Current Ground Speed = ComputeGndSpeedUsingMachAndTrack( Test Speed, Track, 
Altitude k , Wind Speedl, Wind Directionl ) 

Compute the wind speed and direction at the prior altitude. 

InterpolateWindWptAltitude( Wind Profile k ./, Altitudek, Wind Speedl, Wind Directionl ) 
The ground speed at the prior altitude and speed is, 

Prior Ground Speed = ComputeGndSpeedUsingMachAndTrack( Prior Speed, 
GndTrackk-i, Altitudek-i, Wind Speedl, Wind Directionl ) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2 

The distance estimate, dx, is Average Ground Speed * t. 

dx = Average Ground Speed * t / 3600 

Compute the distance required to meet the speed using the previous estimate distance dx. 
Begin by computing the altitude, AltD, at distance dx. 
if (Altitudek > Altitudek-i) AltD = Altitudek 
else 

AltD = (6076 * dx) * tangenU Crossing Angle./,) + Altitudek 
if (AltD >Altitudek-i) AltD = Altitudek 
Compute the average Mach rate. 

MRatel = CasToMach( Crossing Rate cc , Altitudek) 

MRate2 = CasToMach( Crossing Rate cc , AltD ) 
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Test Rate = (MRatel + MRate2) / 2 
t = (Prior Speed - Test Speed) / Test Rate 
The new distance x is DTG k + dx. 

Compute the winds at AltD and distance x. 

InterpolateWindAtDistance( AltD, x, Wind Speed2, Wind Direction2 ) 

The track angle at this point, with GetTrajGndTrk defined in this section, is: 

Track2 = GetTrajGndTrk( x ) 

The ground speed at altitude AltD is then, 

Prior Ground Speed = ComputeGndSpeedUsingMachAndTrack( Prior Speed, Track2, 
AltD, Wind Speed2, Wind Direction2 ) 

Average Ground Speed = (Prior Ground Speed + Current Ground Speed) / 2 

dx = Average Ground Speed * t / 3600 

If there is a TCP prior to dx, compute and insert its speed. 

If the distance is very close to the waypoint, just set the speed. 
if ( (DTG k _i < (DTG k + dx + some small value) ) 
if (\DTG k _j - DTG k - dx \ < some small value) 

Mach k .i = Prior Speed 
Mark TCP k as a Mach segment. 

else 

Compute the speed at the waypoint using v 2 = v 0 2 + 2ax to get v. 

The headwind at the end point is, 

HeadWind2 = Wind Speed2 * cosine(Wind Direction2 - Ground Tracks) 

dx = DTG k .\ - DTG k 

Compute the average Mach rate. 

MRatel = CasToMach( Crossing Rate cc , Altitude k ) 

MRate2 = CasToMach( Crossing Rate cc , Altitude k -i) 
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Test Rate = (MRatel + MRate2) / 2 

The value of Machk-i is computed using function EstimateNextmach, described in 
this section. 

Mach k . i = EstimateNextMach( Test Speed, Current Ground Speed, Prior Speed, 
Head Wind2, Altitude k , dx, Test Rate ) 

Determine if the constraint is met. 

if ( (k- 1) =pc) 

Was the crossing speed met within 0.002 Mach? If not, set this as an error. 
if (\Machp C - Crossing Machp C \ > 0.002) Mark this as an error condition 
Always set the crossing exactly to the crossing speed. 

Mach pc = Crossing Mach pc 
Set the test speed to the computed speed. 

Test Speed = Mach k .j 

Back up the index counter to the next intermediate TCP. 
k = k- 1 

end of if ( (L)TG k _i < (DTG k + dx + some small value) ) 
else 

The constraint occurs between this TCP and the previous TCP. A new VTCP needs 
to be added at this point. 

The along path distance d where the VTCP is to be inserted is: 
d = DTG k + dx 

Save the ground track value at this distance. 

Saved Ground Track = GetTrajGndTrkf d ) 

Insert a new VTCP at location k in the TCP list. The VTCP is inserted between 
TCP/,./ and TCP* from the original list. The function InsertWaypoint should be 
appropriate for the actual data structure implementation of this function. 

InsertWaypoint( k ) 

Update the data for the new VTCP which is now TCP k . 

WptType k = VTCP 
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if (VSegType k = no type) VSegType k = SPEED 


Turn Type k = no turn 
DTG k = d 

The altitude at this point is computed as follows, recalling that the new waypoint is 
TCP k : 

if (Altitude k +i >Altitude k .i) Altitude k = Altitude k . k 

else Altitude k = (6076 * dx) * tangent(Crossing Angle k+I ) + Altitude k +i 

Mach k = Prior Speed 

Mark TCP k as a Mach segment. 

Add the ground track data which must be computed if the new VTCP occurs within a 
turn. The functions WptlnTurn and ComputeGndTrk are described in this sections. 

if (WptInTurn( k )) Ground Track k = ComputeGndTrk( k, d ) 

else Ground Track k = Saved Ground Track 

Compute and add the wind data at distance d along the path to the data of TCP k . 
GenerateWptWindProjile( d, TCP k ) 

Test Speed = Prior Speed 

Since TCP k , has now been added prior to pc, the current constraint counter cc needs 
to be incremented by 1 to maintain its correct position in the list. 

cc = cc + 1 

end of while k > pc. 

Now go to the next altitude change segment on the profde. 
cc = k 

end of while cc > index number of the first waypoint. 

DeltaAngle 

The DeltaAngle function returns angle a, the difference between Anglel and Angle2. The returned 
value may be negative, i.e., -180 degrees > DeltaAngle > 180 degrees. 

a = Angle2 - Anglel 

Adjust "a " such that 0 >a> 360 
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if (a > 180) a = a - 360 


DoT od Acceleration 

The DoTodAcceleration function handles the special case when there is an acceleration to the descent 
Mach at the top-of-descent. This function is invoked from Add Descent Mach Waypoint, which passes in 
the index number for the TOD waypoint, Todlndex, and the Mach value at the TOD, MachAtTod. The 
function will insert the Mach acceleration point into the waypoint list if a valid acceleration point can be 
found. 

Make an initial estimate of the distance to the new Mach value. The function 
TodAccelerationDistance returns the values Valid , k, and dx. 

TodAccelerationDistance( Todldx, MachAtTod, Mach Descent Mach, Valid , k, dx ) 

if (Valid) 

Add the VTCP for the end of the TOD acceleration. 

d = DTG Todidx - dx 

The original ground track will be needed for the new TCP, so save it. 

OldGroundTrack = GetTrajGndTrk( d ) 

Save the wind data at this distance as a temporary TCP. 

GenerateWptWindProfde( d, TemporaryTcp ) 

The new waypoint is downstream of the current value of k. 
k = k + 1 

InsertWaypoint( k ) 

Note that Wpt k is the newly created waypoint. 

WptType k = VTCP 
TurnType k = no turn 

If the new waypoint is not already marked as a special vertical type, mark it as a top-of-descent 
acceleration point. 

if (VSegType k = NONE) VSegType k = TOD acceleration 
DTG k = d 

Calculate the altitude for the new TCP. 

Altitude k = Altituderodidx - (6076 * dx) * tangent(Crossing Angle k + k ) 

Mach k = Mach Descent Mach 

61 



Mach Crossk = Mach Descent Mach 


MachSegmentk = true 

Set the Crossing Rate to the default value of 0.75. 

Crossing Rate k = 0. 75 
Add the appropriate ground track value. 
if (WptInTurn( k )) Ground Tracks = ComputeGndTrk( k, d ) 
else Ground Tracks = OldGroundTrack 
Copy the wind data from TemporaiyTcp into Wpt k . 
end of if (Valid) 

else set an error for being unable to accelerate to the descent Mach value. 

EstimateN extC as 

EstimateNextCas is an iterative function to estimate the CAS value, CAS, at the next TCP. Note that 
there is no closed-form solution for this calculation of CAS. The input variable names described in this 
function are from the calling routine and are, in order, the target CAS value, Test CAS, the ground speed 
at the estimation starting point, Current Ground Speed, an estimation limiting flag, No Limit Flag, the 
CAS at the estimation starting point, Prior CAS, the head wind at the estimation starting point, Head 
Wind, the altitude at the estimation starting point, Altitude, the distance from the estimation starting point 
to the point where the CAS is to be estimated, Distance, and the deceleration rate to be used in this 
estimation, CAS Rate. Also, the input deceleration value must be greater than 0, CAS Rate > 0. The 
function returns the estimated CAS value. 

Guess CAS = Test CAS 

Set up a condition to get at least one pass. 

d = -10 * Distance 

size = 1.01 * (Prior CAS - Guess CAS) 
count = 0 

if ( (Distance > 0) and (CAS Rate > 0) ) 

Iterate a solution. The counter count is used to terminate the iteration if the distance estimation 
does reach a solution within 0.001 nmi. 

while ( (\Distance - d\ > 0.001) and (count < 10) ) 

if (Distance > d) Guess CAS = Guess CAS - size 

else Guess CAS = Guess CAS + size 
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size = size / 2 


The estimated time t to reach this speed, 
t = (Guess CAS - Test CAS) / CAS Rate 
The new ground speed, 

Gs2 = CasToTas Conversion( guess, Altitude ) - Head Wind 
d = ((Current Ground Speed + Gs2) / 2) * (t / 3600) 
count = count + 1 
end of the while loop 
Limit the computed CAS, if necessary. 

if (( NoLimit = false) and (Guess CAS > Prior CAS)) Guess CAS = Prior CAS 
return Guess CAS 

EstimateNextMach 

EstimateNextMach is an iterative function to estimate the Mach value, Mach, at the next TCP. Note 
that there is no closed-form solution for this calculation. The input variable names described in this 
function are from the calling routine. Also, the input deceleration value must be greater than 0, Mach 
Rate > 0. 

Mach = Test Speed 

Set up a condition to get at least one pass. 
d = -10 * dx 

size = 1.01 * (Prior Speed - Test Speed) 
count = 0 

if ( (dx > 0) and (Test Rate > 0)) 

Iterate a solution. The counter count is used to terminate the iteration if the distance estimation 
does reach a solution within 0.001 n.mi. 

while ( ( d - dx \ > 0.001) and (count < 10)) 

if (d > dx) Mach = Mach - size 

else Mach = Mach + size 

size = size / 2 

The estimated time t to reach this speed. 
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t = (Mach - Test Speed) / Test Rate 
The new ground speed, 

CAS = MachToCas( Mach, Altitude ) 

Gs2 = CasToTas Conversion( CAS, Altitude ) - Head Wind2 
d = ((Current Ground Speed + Gs2) / 2) * (t / 3600) 
count = count + 1 
end of the while loop 
Limit the computed Mach, if necessary. 
if (Mach > Prior Speed) Mach = Prior Speed 
GenerateWptWindProfile 

The function GenerateWptWindProfile is used to compute new wind profile data. This function is a 
double-linear interpolation using the wind data from the two bounding input waypoints to compute the 
wind profile for a new VTCP, TCP *. The interpolations are between the wind altitudes from the input data 
and the ratio of the distance d at a point between TCP,./ and TCP, and the distance between TCP,./ and 
TCP, E.g., 

- Find the two bounding input waypoints, TCP,] and TCP „ between which d lies, e.g., 
TCP,, >d> TCP, 

- Using the altitudes from the wind profile of TCP „ compute and temporarily save the wind data at 
these altitudes using the wind data from TCP,, (e.g., Wind Speed Temporary, Aititudei)- 

- Compute the wind speed and wind direction for each altitude using the ratio r of the distances. 
Assuming that the difference between DTG,./ and DTG, f 0, and that DTG„ > DTG, 

r = (DTG,., - d) / (DTG,, - DTG,) 

Iterate the following for each altitude in the profile. 

Wind Speedy Aititudei = ( (1 - r) * Wind Speed Temporary, Aititudei) + (r * Wind Speed, Aititudei) 
a = DeltaAngle(Wind Direction Temporary, Aititudei, Wind Direction, Aititudei) 

Wind Direction k , Aititudei = Wind Direction k , Aititudei + O' * a) 
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Figure 7 is an example of the computation data for the wind computation at a 9,000 ft altitude. In this 
example, TCP has wind data at 10,000 and 8,000 ft and TCP, has wind data at 9,000 ft. 


Wind Data^j, -\- 


10000ft 


Wind Data} j 4-/ 
8000ft ’ 


\ (1 - r) * Wind Data^j 

^ computed for 9000 ft 

H “ H 

DTG/_] -d 


r * Wind Data; 


DTG i _ 1 -DTG t 


Wind Data 
9000ft 


l, 


TCPj_j TCP k TCPj 

Figure 7. Example for computing a single wind data altitude. 


GetT raj ectoryData 

The GetTrajecton’Data function computes the trajectory data at the along-path distance equal to d and 
saves these data in a temporary TCP record. The function uses a linear interpolation based on the DTG 
values of the two TCPs bounding this distance and the distance d to compute the trajectory data at this 
point. 

GetTrajGndTrk 

The GetTrajGndTrk function computes the ground track at the along-path distance, distance, 
if ( (distance < 0) or (distance > DTGf irstwavp omt) ) Ground Track = Ground Trackf, rstwayp0 int 

else 


Find where distance is on the path. 

i = index number of the last waypoint 

while (distance > DTG j i = i -1 

if (distance = DTG,) Ground Track = Ground Track , 

else 


x = DTGi - DTG i+1 
if (x <0) r = 0 

else r = (distance - DTG i+ i) / x 
if (r > 1) r = 1 

dx = (1 - r) * DeltaAngle( Ground Track „ Ground Track i+I ) 
Ground Track = Ground Trackj + dx 
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HandleDescentAccelDecel 

The function HandleDescentAccelDecel is designed to handle the special case of a Mach acceleration 
in the descent where the first CAS crossing restriction cannot be met. The calling program provides as 
input and retains the subsequent outputs for the following variables: Caslndex, CruiseMach, 
MachCasModified, DescentMach, and MachCas. The variable Caslndex is the index value in the TCP list 
for the first CAS constrained waypoint. The variable CruiseMach is the last Mach crossing restriction 
value prior to the first CAS segment. The variable MachCasModified is a flag returned by this function if 
the DescentMach or MachCas values are changed. The variables DescentMach and MachCas are the 
planned descent Mach and planned Mach-to-CAS transition CAS, respectively, and these values may be 
modified by this function. 

Initialize variables. 

i = 0 

z = 0 

find = false 

MachCasModified = false 

Perform up to two iterations to calculate any required Mach or CAS change in the descent. 

while ( (fini = false) and (i < 2)) 

Calculate z at the descent Mach and the Mach-to-CAS CAS. 
z = MachCasTransitionAltitude( MachCas, DescentMach ) 

Determine if z is below the CAS crossing restriction. 


if (z < A l titllde Caslndex) 


Set the CAS to the value at this altitude, knowing the crossing restriction can't be met. 
MachCas = MachToCas( DescentMach, Altitudecasindex) 
else if ( z> Altitude Cross fi rstwa y poin ,) 

Set the Mach to the descent CAS at the cruise altitude. 
m = CasToMach( MachCas, Altitude first waypoint) 
if (m > CruiseMach) DescentMach = m 
if ( MachCas <Crossing CAScasinda ) 

MachCas = Crossing CAScasindex 
i = i + 1 
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else fini = true 

end of while ( (fini = false) and (i < 2)) 

Find the TOD TCP. 
fini = false 
Todlndex = 0 

i = index number of the first waypoint 

while ((i < index number of the last waypoint) and (fini = false)) 
if ( ( Altitudei < Altitude f irs , waypoint) or (Crossing CAS , > 0) ) 
if (( Altitudei f Altitudefirst wa ypoini)) Todlndex = i - 1 
else Todlndex = i 
fini = true 
i = i + 1 

end of while (fi < index number of the last waypoint) and (fini = false)) 

Calculate the entire decent distance. 

d ~ TTTGfodlndex ~ DTGcasIndex 

Estimate the distance, Daccel, to the new Mach value. 

TodAccelerationDistance( Todlndex, CruiseMach, MachDescentMach, Valid, Accellndex, Daccel ) 
Estimate the distance, Ddecel, to the CAS crossing speed. 

BodDecelerationDistance( Caslndex, z, Mach Transition CAS, Ddecel ) 
fini = false 
m = DescentMach 

The nominal speed values won't work, there is insufficient distance to obtain the acceleration and then 
slow to the crossing speed. Iterate until a solution is found. 

while ((fini = false) and (d < (Daccel + Ddecel))) 

Iterate the solution. 

Slightly change the Mach and then find the CAS. 

7 7i = m - 0. 002 
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if (m < Cruise Mach) 


m = Cruise Mach 
fini = true 

Estimate the distance to the new Mach value. 

TodAccelerationDistance( Todlndex, Cruise Mach, m, Valid, Accellndex, Daccel ) 

Find the altitude where the acceleration ends. 

z = Crossing Altitude j\ rs t wavp omt - (Daccel / d) * (Crossing Altitude first waypoint - 
Crossing Altitude Cas index) 

CAS = MachToCas( m, z ) 

Estimate the distance to the CAS crossing speed. 

BodDecelerationDistance( Caslndex, z, CAS, Ddecel ) 

if (d> (Daccel + Ddecel)) 

fini = true 

Modify the descent Mach and CAS values, 
modified = true 
DescentMach = m 

Add a buffer to the CAS so that subsequent Mach-to-CAS calculation won't cause an error. 
MachCas = CAS + 0.1 
end of if (d > (Daccel + Ddecel)) 

InterpolateWindAtDistance 

The function InterpolateWindAtDistance is used to compute the wind speed and direction at an 
altitude, Altitude, for a specific distance, Distance, along the path. This function is a linear interpolation 
using the wind data from the input waypoints that bound the along-path distance. 

Find the bounding input waypoints. 

iO = index number of the first waypoint 

j = index number of the first waypoint 

fini = false 

if (Distance < 0) Distance = 0 
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while ( (fini = false) and (j < (index number of the last waypoint - 1) ) ) 
if ( (TCPj = input waypoint) and (DTGj > Distance) ) iO =j 
if(DTGj < Distance ) fini = true 
end of the while loop 
il = iO + 1 

j = H 

fini = false 

while ( (fini = false) and (j < index number of the last waypoint) ) 
if ( (TCPj = input waypoint) and (DTGj < Distance) ) 

il =j 

fini = true 
end of if 

j=j + l 

end of the while loop 

if (il > index number of the last waypoint) il = index number of the last waypoint 

if (iO = il) InterpolateWindWptAltitude( TCP i0 , Altitude ) 

else 

Interpolate the winds at each waypoint. 

InterpolateWindWptAltitude( TCP i0 , Altitude ), returning SpdO and DirO 
InterpolateWindWptAltitude( TCPu, Altitude ), returning Spdl and Dirl 
Interpolate the winds between the two waypoints. 
r = (DTGjo - Distance) / (DTG,o - DTGfi 
Wind Speed = ( (1 - r) * SpdO) + (r * Spdl) 
a = DeltaAngle( DirO, Dirl ) 

Wind Direction = DirO + (r * a) 
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InterpolateWindWptAltitude 

The function InterpolateWindWptAltitude is used to compute the wind speed and direction at an 
altitude, Altitude, for a specific TCP. This function is a linear interpolation using the wind data from the 
current TPC. 

Find the index numbers, pO and pi, for the bounding altitudes. 
pO = 0 
pi =0 

for (k = 1; k< Number of Wind Altitudes k = k + 1) 
if (Wind Altitude ii k < Altitude) pO = k 
if ( (Wind Altitude k > Altitude) and (pi = 0) ) pi = k 
if (pi = 0) pi = Number of Wind Altitudes 

If Altitude = Wind Altitude p o or if Altitude = Wind Altitude p i then the wind data from that point is 
used. Otherwise, Altitude is not at an altitude on the wind profile of TCP,, i.e., z = Wind Altitude/, k , 
then: 

if (Wind Altitude p i < Wind Altitude p o) r = 0 

else r = (Altitude - Wind Altitude p0 ) / (Wind Altitude pI - Wind Altitude p0 ) 

Wind Speed = ( (1 - r) * Wind Speed p0 ) + (r * Wind Speed p i) 
a = DeltaAngle(Wind Direction p o, Wind Direction p i) 

Wind Direction = Wind Direction p0 + (r * a) 

MachCasT ransitionAltitude 

The function MachCasTransitionAltitude is used to compute the altitude where the input Mach, Mach, 
and CAS, Cas, values would be equivalent 

z = (1 - (((((0.2 * ((Cas/661.48) 2 ) + if 5 ) - 1) / (((0.2 * (Mach 2 ) + if 5 ) - l)f 19026 )) / 0.00000687535 
return the value of z. 

RadialRadiallntercept 

The function RadialRadiallntercept determines if two place-and-radial sets, each defined by a latitude, 
a longitude, and a track angle, will intersect and if so, calculates the latitude and longitude of the intercept 
point. Inputs are values of latitude, Latitude, longitude, Longitude, and angle, Angle', one set of each for 
the two place-and-radial sets. If a valid intercept can be calculated, then the intercept point's latitude and 
longitude are output, NewLatitude and NewLongitude, and the function returns a valid indication. 
Otherwise, the function returns an invalid indication. 

Calculate the distance and the track angle between the two input positions. 
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distance i, 2 = arccosine( sine(Latitude 1 ) * sine (Latitude 2 ) + cosine(Latitudei) * cosinet Latitude 2 ) * 
cosine(Longitudei - Longitude 2 ) ) 

track 12 = arctangent2( sine(Longitude 2 - Longitude 1 ) * cos in e ( La titude 2 ) , cosine(Latitudei) * 

sine(Latitude 2 ) - sine(Latitude 1 ) * cosine(Latitude 2 ) * cosine(Longitude 2 - Longitude!)) 

Check for error in the intercept calculation. 

error = false 

track 1 = Angle! - track i, 2 + 90 

Adjust tracki such that 0 > track 1 > 360 

track 2 = Angle 2 - tracki, 2 + 90 

Adjust track 2 such that 0 > track 2 > 360 

Determine the quadrant. 

angi = track 2 + 180 

Adjust angi such that 0 > angi > 360 

if ( ( \DeltaAngle(trackl , track2)\ < 2) or ( \DeltaAngle(trackl , angi) \ < 2)) 

Determine if the angles are really 180 degrees apart. 

ang 2 = Angle 2 + 180 

Adjust ang 2 such that 0 > ang 2 > 360 

ang } = Del taAngleLAngle /, ang 2 ) 

ang 4 = Del taA ngle (Angi e /, track j 2 ) 

if ((\ang3\ > 2) or (\ang4\ > 2)) error = true 

if (error = false) 

RelativeLatLong(Latitude h Longitude 1 , track i 2 , distance i, 2 / 2, NewLatitude, NewLongitude) 

else 

Determine the quadrant. 
if (tracki < 90) quadrant 1 = 1 
else if (tracki < 180) quadrantl = 2 
else if (tracki < 270) quadrantl = 3 


71 



else quadrantl = 4 
if (track 2 < 90) quadrant2 = 1 
else if (track 2 < 180) quadrantl = 2 
else if (track 2 < 270) quadrantl = 3 
else quadrant 2 = 4 
if (quadrantl = 1) 

if ((quadrant2 = 2) or (quadrant2 = 3 )) error = true 
if (( quadrant 2 = 1) and (chktkl < chktk2)) error = true 
else if (quadrantl = 2) 

if (( quadrant2 = 1) or (quadrant2 = 4)) error = true 
if ((quadrant2 = 2) and ( chktkl > chktk2)) error = true 
else if (quadrantl = 3) 

if ((quadrant2 = 1) or (quadrant2 = 2) or ( quadrant 2 = 4)) error = true 
if (track: > track 2 ) error = true 

else 

if (( quadrant 2 = 1) or (quadrant2 = 2) or (quadrantl = 3f) error = true 
if (track: < track 2 ) error = true 
if (error = false) 

trx : = \Angle: - track] 2 \ 

Adjust trx : such that 0 > trx: > 360 
trx 2 = \Angle 2 - (track: _ 2 + 180) \ 

Adjust trx 2 such that 0 > trx 2 >360 
if (trx: > 180) trx: = 360 - trx : 
if (trx 2 > 180) trx 2 = 360 - trx 2 
ang 5 = 180 - trx: - trx 2 

if ((angs = 0) or ((ang 5 -180) = 0) or (distance:? = 0 )) error = true 
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if (error = false) 


distance 2 = distance u * sine( trx 2 ) / sine(ang 5 ) 
if (distance 2 < 0) distance 2 = - distance 2 
if (distance 2 > max intercept range) error = true 

else RelativeLatLong(Latitude i, Longitude i, Angle i, distance 2 , NewLatitude, 
NewLongitude) 

if (error) return false 

else return true 

RelativeLatLon 

The function RelativeLatLon computes the latitude and longitude from input values of latitude, 
BaseLat, longitude, BaseLon, angle, Angle, and range, Range. 

if (Angle = 180) Latitude = -Range / 60 + BaseLat 

else Latitude = ( (Range * cos(Angle) ) / 60) + BaseLat 

if ( (BaseLat = 0) or (BaseLat = 180) ) Longitude = BaseLon 

else if (Angle = 90) Longitude = BaseLon + Range / (60 * cos (BaseLat) ) 

else if (Angle = 270) Longitude = BaseLon - Range / (60 * cos (BaseLat) ) 

else 


rl = tangent(45 + 0.5 * Latitude) 
r2 = tangent(45 + 0.5* BaseLat) 

if ( (rl = 0) or (r2 = 0) ) Longitude = 20, just some number, mark this as an error condition, 
else Longitude = BaseLon + (180/pi * (tangent (Angle)* (log(rl) - log(r2)))) 

T odAccelerationDistance 

The TodAccelerationDistance function estimates the distance required for the special case of an 
acceleration from the top-of-descent Mach to the descent Mach at the top-of-descent. This function is 
invoked from HandleDescentAccelDecel and DoTod Acceleration, which passes in the index number for 
the TOD waypoint, Todlndex, and the Mach value at the TOD, MachAtTod. The function returns a 
validity flag to indicate if a TOD acceleration is valid. Valid, and if valid, the indices in the TCP list 
where the acceleration occurs, Accellndex, and the distance from the index point of the acceleration, 
Distance. 

Perform an initialization of flags and counters. 
fini = false 
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skip = true 
k = Todlndex 

Make an initial guess of the distance to the new Mach value. 

Descent Speed = Mach Descent Mach 

Mach Ratei = CasToMach( 0. 75 kt / sec, Altitude Todindex) 

Compute the time required to do the deceleration. 
t = (Mach Descent Mach — MachAtTod ) / Mach Ratei 
Compute the wind speed and direction at the current altitude. 

InterpolateWindWptAltitude( Wind Profile To dindex> Altitude Todindex, Wind Speed, Wind Direction ) 

Get the ground track at the current point. 

if ( WptInTurn( Waypoint To dindex) ) track = Ground Trackfodindex + 1 
else track = Ground Track To dindex 

TOD Ground Speed = ComputeGndSpeedUsingMachAndTrackf MachAtTod, track, Altitude To dindex, 
Wind Speed, Wind Direction ) 

Descent Ground Speed = ComputeGndSpeedUsingMachAndTrack( Mach Descent Mach, track, 
Altitude To dindex, Wind Speed, Wind Direction ) 

The average ground speed is as follows: 

Average Ground Speed = ( TOD Ground Speed + Descent Ground Speed ) / 2 
The distance estimate, dx, is Average Ground Speed * t with a conversion to nm. 
dx = Average Ground Speed * t / 3600 

Now compute better estimates, doing this twice to refine the estimation. 
for ( i = 1; i <2; i = i + 1 ) 
skip = false 

Determine if this distance is beyond the next downstream waypoint. 
k = Todindex 
d DTGfodlndex ~ dx 

while ( ( k < ( index number of the last waypoint - 1 ) ) and ( DTGk+i > d ) ) 
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if ( ( kf Todlndex ) and ( Crossing Rate * > 0 ) ) skip = True; 

k = k + 1 

Compute the wind speed and direction at the new altitude. 

InterpolateWindWptAltitude( Waypoint k , Altitude k , Wind Speed, Wind Direction ) 

The ground speed at the this point is: 

Descent Ground Speed = ComputeGndSpeedUsingMachAndTrack( Mach Descent Mach, Ground 
Track h Altitude^ Wind Speed, Wind Direction ) 

The average ground speed is: 

Average Ground Speed = ( TOD Ground Speed + Descent Ground Speed ) / 2 
The distance, dx, is: 
dx = Average Ground Speed * t / 3600 
If there is a valid deceleration point, add it. 

Valid = not skip 
Accellndex = k 
Distance = dx 
WptlnTurn 

The WptlnTurn function simply determines if the waypoint is between a turn-entry TCP and a turn-exit 
TCP. If this is true, then the function returns a value of true, otherwise it returns a value of false. 

fini = false 

within = false 

j = i + 1 

while ((fini = false) and (j < (index number of the last waypoint))) 
if (TurnTypej = turn-entry) fini = true 
else if (TurnTypej = turn-exit) 
fini = true 
within = true 

j=j + l 

return within 
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Summary 

The algorithm described in this document takes as input a list of waypoints, their trajectory-specific 
data, and associated wind profile data. A full 4D trajectory can then be generated by the techniques 
described. A software prototype has been developed from this documentation. An example of the data 
input and the prototype-generated output is provided in the Appendix. 
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Appendix Example Data Sets 

Input Trajectory Data 

An example input trajectory data set is provided in Table Al. 

The descent Mach is 0.82. The Mach-to-CAS transition speed for this example is 310 knots. Note that 
Waypoint-18 is the runway threshold at a 50 ft crossing height. No RF turns are shown. 


Table Al. Example of trajectory input data. 


Identifier 

Latitude 

Longitude 

Crossing 

Altitude 

Crossing 

Angle 

Crossing 

CAS 

Crossing 

Mach 

Crossing 

Rate 

Waypoint-01 

31.87476 

-103.244 

37000 

0 

0 

0.78 

0 

Waypoint-02 

32.48133 

-99.8635 

0 

0 

0 

0 

0 

Waypoint-03 

32.20548 

-98.9531 

0 

0 

0 

0 

0 

Waypoint-04 

32.19398 

-98.6621 

0 

0 

0 

0 

0 

Waypoint-05 

32.17042 

-98.113 

0 

0 

0 

0 

0 

Waypoint-06 

32.15959 

-97.8777 

0 

0 

0 

0 

0 

Waypoint-07 

32.34026 

-97.6623 

0 

0 

0 

0 

0 

Waypoint-08 

32.46908 

-97.5079 

0 

0 

0 

0 

0 

Waypoint-09 

32.64444 

-97.2967 

11700 

3.0 

0 

0 

0 

Waypoint-10 

32.71448 

-97.2119 

11000 

1.1 

0 

0 

0 

Waypoint- 1 1 

32.74948 

-97.1695 

0 

0 

0 

0 

0 

Waypoint- 12 

32.97496 

-97.1783 

0 

0 

0 

0 

0 

Waypoint- 13 

33.10724 

-97.1754 

5300 

2.3 

220 

0 

0.5 

Waypoint- 14 

33.10658 

-97.0537 

4300 

1.8 

0 

0 

0 

Waypoint- 15 

33.03645 

-97.0541 

0 

0 

0 

0 

0 

Waypoint- 16 

33.00561 

-97.0542 

2400 

3.1 

170 

0 

0.5 

Waypoint- 17 

32.95953 

-97.0544 

1495 

3.0 

127 

0 

0.75 

Waypoint- 18 

32.91582 

-97.0546 

660 

3.0 

127 

0 

0.75 
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Input Wind Data 

An example wind speed data set is provided in Table A2. 


Table A2. Example of wind speed input data. 




Wind 

Wind 

Identifier 

Altitude 

Speed 

Direction 

Waypoint-01 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-02 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-03 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-04 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-05 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-06 

0 

20 

180 


10000 

50 

270 


20000 

60 

340 


40000 

70 

350 

Waypoint-07 

0 

20 

160 


10000 

50 

240 


20000 

60 

320 


40000 

70 

330 
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Table A2 (continued). Example of wind speed input data. 


Identifier 

Altitude 

Wind 

Speed 

Wind 

Direction 

Waypoint-08 

0 

20 

160 


10000 

50 

240 


20000 

60 

330 


40000 

70 

340 

Waypoint-09 

0 

20 

160 


10000 

50 

240 


20000 

60 

330 


40000 

70 

340 

Waypoint- 10 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 1 1 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint-12 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 13 

0 

20 

160 


10000 

50 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 14 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 
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Table A2 (continued). Example of wind speed input data. 


Identifier 

Altitude 

Wind 

Speed 

Wind 

Direction 

Waypoint- 15 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 16 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 17 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 

Waypoint- 18 

0 

20 

160 


10000 

40 

240 


20000 

50 

330 


40000 

60 

340 


Output Trajectory Data 

An example of the data available from this trajectory algorithm is provided in Table A3. Not shown, 
but also available, are the latitude and longitude data for each TCP. 


Table A3. Example of the trajectory output data. 


TCP type 

Identifier 

Altitude 

Mach 

CAS 

Mach 

Segment 

Ground 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint-01 

37000 

0.78 

252.5 

true 

450.7 

77.1 

366.06 

3214.8 

Turn-entry 


37000 

0.78 

252.5 

true 

450.7 

77.1 

192.89 

1831.4 

Input 

Waypoint-02 

37000 

0.78 

252.5 

true 

469.9 

93.3 

190.64 

1813.8 

Turn-exit 


37000 

0.78 

252.5 

true 

487.5 

109.5 

188.39 

1796.9 

Turn-entry 


37000 

0.78 

252.5 

true 

487.5 

109.5 

142.90 

1461.0 

Input 

Waypoint-03 

37000 

0.78 

252.5 

true 

478.6 

101 

141.68 

1451.9 

Turn-exit 


37000 

0.78 

252.5 

true 

469.1 

92.6 

140.46 

1442.6 
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Table A3 (continued). Example of the trajectory output data. 


TCP type 

Identifier 

Altitude 

Mach 

CAS 

Mach 

Segment 

Ground 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint-04 

37000 

0.78 

252.5 

true 

469.1 

92.8 

126.90 

1338.6 

VTCP 


37000 

0.78 

252.5 

true 

469.3 

93 

125.46 

1327.5 

VTCP 


36306 

0.82 

271.2 

true 

494.5 

93 

123.28 

1311.2 

VTCP 


30337 

0.82 

310 

false 

509.6 

93 

104.53 

1176.8 

Input 

Waypoint-05 

28569 

0.793 

310 

false 

497.2 

93 

98.98 

1137.1 

Turn-entry 


25777 

0.751 

310 

false 

478.5 

93 

90.21 

1072.4 

Input 

Waypoint-06 

24818 

0.737 

310 

false 

446.6 

69.1 

87.20 

1048.9 

Turn-exit 


23858 

0.723 

310 

false 

415.4 

45.2 

84.19 

1023.8 

Input 

Waypoint-07 

19976 

0.672 

310 

false 

393.4 

45.3 

72.00 

915.2 

Input 

Waypoint-08 

16474 

0.629 

310 

false 

404.6 

45.4 

61.00 

816.0 

Input 

Waypoint-09 

11700 

0.576 

310 

false 

409.4 

45.5 

46.01 

683.4 

VTCP 


11432 

0.574 

310 

false 

408.5 

45.5 

43.71 

663.1 

Input 

Waypoint- 10 

11000 

0.524 

284.6 

false 

378.1 

45.5 

40.01 

629.3 

VTCP 


11000 

0.519 

282 

false 

375.1 

45.5 

39.65 

625.8 

Turn-entry 


10811 

0.507 

276.4 

false 

368.4 

45.5 

38.87 

618.3 

Input 

Waypoint- 1 1 

10382 

0.479 

262.9 

false 

340.6 

21.8 

37.12 

600.5 

VTCP 


10000 

0.453 

250 

false 

324.7 

19.3 

35.55 

583.5 

Turn-exit 


9954 

0.452 

250 

false 

308.9 

358.1 

35.36 

581.4 

Input 

Waypoint- 12 

7105 

0.429 

250 

false 

307.7 

1.1 

23.69 

445.1 

VTCP 


6474 

0.424 

250 

false 

307.3 

1.1 

21.10 

414.8 

Turn-entry 


5793 

0.391 

233.1 

false 

286.5 

1.1 

18.31 

381.0 

Input 

Waypoint- 13 

5300 

0.366 

220 

false 

270 

45.7 

16.29 

354.9 

Turn-exit 


4909 

0.363 

220 

false 

245 

90.3 

14.27 

326.6 

Turn-entry 


4556 

0.361 

220 

false 

242 

90.3 

12.42 

299.3 

Input 

Waypoint- 14 

4300 

0.359 

220 

false 

215.4 

135.3 

11.08 

278.2 

VTCP 


3987 

0.357 

220 

false 

204.1 

164.4 

10.21 

263.2 

Turn-exit 


3831 

0.35 

215.9 

false 

197 

180.3 

9.74 

254.7 

Input 

Waypoint- 15 

3009 

0.305 

191.2 

false 

170.7 

180.2 

7.24 

205.8 

Input 

Waypoint- 16 

2400 

0.268 

170 

false 

148.8 

180.2 

5.39 

164.1 

VTCP 


2140 

0.267 

170 

false 

148.9 

180.2 

4.65 

146.2 


Table A3 (continued). Example of the trajectory output data. 
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Mach 

Ground 




TCP type 

Identifier 

Altitude 

Mach 

CAS 

Segment 

Speed 

Track 

DTG 

TTG 

Input 

Waypoint- 17 

1495 

0.197 

127 

false 

105.5 

180.2 

2.62 

88.9 

Input 

Waypoint- 18 

660 

0.194 

127 

false 

106.9 

180.2 

0.00 

0.0 
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